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ABSTRACT 


Middle  ultraviolet  spectra  from  1800  to  340oA  are 
analyzed.  These  spectra  were  obtained  from  the  March  1990 
rocket  flight  of  the  NPS  MUSTANG  instrument  over  the  altitudes 
105km  to  315km.  The  data  were  compared  with  computer 
generated  synthetic  spectra.  A  least  sguares  fitting 
procedure  was  developed  for  this  purpose.  Each  data  point  was 
weighted  using  the  standard  deviation  of  the  means.  Synthetic 
spectra  were  generated  for  the  following  emissions:  N2 
Vegard-Kaplan;  Nj  Lyman-Birge-Hopf ield;  NO  7,  S,  and  e;  01 
2972A,  Oil  247oA;  and  Nil  2143A.  Altitude  profiles  for  the 
emissions  were  obtained.  Tentative  identification  was  made  of 
the  OIII  2853A  emission.  A  comparison  of  VK  and  LBH  profiles 
demonstrates  the  process  of  N2  A-state  quenching  by  atomic 
oxygen . 
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I. 


INTRODUCTION 


The  composition  of  the  earth's  atmosphere  is  important 
to  many  areas  of  science,  such  as  weather  prediction, 
satellite  operations,  and  communications.  Properties  of  the 
upper  atmospheric  region  known  as  the  ionosphere  are  of 
specific  interest  to  the  military,  as  they  strongly  influence 
such  systems  as  high  frequency  (HF)  radio  communications, 
over-the-horizon  radar  (OTH) ,  ballistic  missile  early  warning 
(SDI),  and  the  Ground  Wave  Emergency  Network  (GWEN).  The 
electron  density  profile,  a  plot  of  electron  density  versus 
altitude  for  a  given  time  and  location,  is  needed  for 
effective  use  of  the  ionosphere.  Electron  density  is  a  fav^tor 
that  varies  not  only  with  altitude,  but  with  time  and 
geographical  location  as  well.  In  1936,  the  Joint  Chiefs  of 
Staff  prioritized  the  determination  of  ionospheric  electron 
density  as  fifth  on  a  list  of  43  critical  geophysical 
parameters  for  inclusion  in  a  defense  satellite  (The  Joint 
Chiefs  of  Staff, 1986). 

Currently,  electron  density  profiles  are  determined  using 
ground-based  ionosonde  stations.  This  is  very  effective 
locally,  but  is  inherently  inaccurate  for  global  forecasting 
since  there  are  only  about  20  stations  worldwide.  A 
satellite-based  method  of  forecasting  the  ionosphere  would  be 
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globally  effective,  but  to  place  an  ionosonde  in  orbit  is 
unfeasible  due  to  size  and  power  limitations.  One  alternative 
to  ionosondes  is  to  infer  electron  density  profiles  from 
spectrographic  measurements  of  the  atmospheric  airglow.  An 
operational  spectrograph  is  small  and  relatively  low-powered; 
it  can  easily  be  included  in  a  defense  satellite  package.  The 
Naval  Research  Lab  and  the  Naval  Postgraduate  School  have  a 
joint  spectrographic  project  scheduled  to  be  included  in  a 
satellite  to  be  launched  in  1995.  The  NRL  project  is  a  high 
resolution  airglow  and  aurora  spectrograph  (HIRAAS)  and  the 
accompanying  NPS  effort  is  the  middle  ultraviolet  spectrograph 
(MUSTANG) .  The  inclusion  of  the  duo  represents  a  significant 
improvement  in  the  global  determination  of  electron  density 
profiles  over  ground-based  ionosondes. 

The  inference  of  electron  densities  from  emission 
observations  by  the  HIRAAS/MUSTANG  is  performed  by  relating 
measurements  of  natural  ultraviolet  emissions  to  the  neutral 
species'  densities.  A  rocket-launched  HIRAAS/MUSTANG 
experiment  was  first  tested  in  March  1990.  The  data  from  the 
launch  is  the  topic  of  this  thesis. 

A.  THESIS  OBJECTIVES 

This  thesis  focuses  on  the  analysis  of  data  obtained  from 
the  March  1990  flight.  The  analysis  extends  methods  used 
previously  on  subsets  of  the  data  to  the  entire  observed 
range,  ISOOA  to  340oA.  All  known  synthetic  models  with 
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emissions  in  this  range  are  used  in  the  fitting.  The  result 
of  the  fitting  algorithm  is  a  set  of  relative  scale  factors 
for  the  various  spectra  by  altitude.  The  relative  scale 
factors  are  converted  into  absolute  terms  and  plotted  against 
altitude.  These  absolute  plots  of  intensity  and  column 
density  can  be  used  to  infer  an  electron  density  profile  (see 
e.g.,  Meier, 1991;  McCoy, 1985).  A  secondary  goal  is  to 
determine  the  feasibility  of  an  approximation  for  the  N2 
Second  Positive  spectrum  given  the  unavailability  of  a 
synthetic  spectrum.  The  determination  of  the  Second  Positive 
spectrum  would  allow  an  estimation  of  photoelectron  flux, 
since  it  is  produced  by  transitions  from  a  state  which  can 
only  be  populated  by  collisions. 

B.  THESIS  OUTLINE 

The  thesis  is  divided  into  six  chapters.  Chapter  II  gives 
general  background  information  on  the  descriptions  and 
designations  of  the  atmosphere  and  existing  emissions.  A 
basic  explanation  of  atomic  and  molecular  spectra  is  also 
summarized  in  this  chapter.  Chapter  III  addresses  the 
specific  emission  features  in  the  wavelength  range  of 
interest.  Chapter  IV  presents  a  summary  of  the  experiment  and 
the  subsequent  data  reorganization.  Chapter  V  gives  a 
complete  treatment  of  analysis  of  the  data  and  results  of  the 
study.  Chapter  VI  concludes  the  paper  with  a  summary  of 
findings  and  suggestions  for  future  study.  The  Appendices 
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contain  the  computer  programs  which  generated  the  synthetic 
spectra  and  the  fitting  programs,  as  well  as  plots  of  the 
final  fits  to  the  data. 
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II.  BACKGROUND 


A.  THE  ATMOSPHERE 

1.  General  Description 

The  Earth's  atmosphere  is  designated  by  horizontal  layers, 
and  these  divisions  are  determined  by  one  of  two  primary 
methods.  One  method  designates  the  layers  based  upon 
temperature  characteristics;  the  other  uses  density  to 
categorize  the  strata.  The  layers  and  a  typical  temperature 
profile  are  shown  in  Figure  2-1. 
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Figure  2-1  Temperature  profile  and  corresponding  atmospheric 
layers . 
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With  respect  to  the  temperature  designations,  the 
troposphere  is  the  region  closest  to  the  Earth's  surface  where 
temperatures  decrease  with  altitude,  typically  from  about  300K 
at  the  surface  to  220K  at  10  km.  The  point  at  which 
temperature  begins  to  rise  again  is  called  the  tropopause  and 
it  indicates  the  start  of  the  next  region,  the  stratosphere . 
The  stratosphere  is  a  relatively  warm  layer  due  to  absorption 
of  solar  radiation (primarily  2000-3000A)  by  ozone.  The  top  of 
the  stratosphere  is  marked  by  the  stratopause.  At  this  point 
the  mesosphere  begins.  The  temperature  again  begins  to 
decrease  with  altitude  until  a  low  of  about  lOOK  occurs  at 
around  85  kilometers.  This  is  the  mesopause ,  which  marks  the 
beginning  of  the  next  layer,  the  thermosphere .  In  the 
thermosphere,  temperature  rises  steadily  with  increasing 
altitude  until  reaching  an  equilibrium  temperature  between 
lOOOK  and  2000K,  depending  on  solar  activity. 

Studying  the  atmosphere  in  terms  of  density  gives  rise  to 
a  different  set  of  atmospheric  layers.  Under  this  convention, 
the  lowest  80  kilometer  layer  of  the  atmosphere  is  called  the 
homosphere .  In  this  region,  vertical  mixing  of  the  air 
currents  provides  a  uniform  relative  density  ratio  between 
various  species.  The  primary  components  are  nitrogen,  oxygen, 
and  argon  in  a  ratio  of  78:21:0.9.  The  remaining  0.1  percent 
is  primarily  helium,  hydrogen,  and  carbon  dioxide.  The  20  km 
from  80  to  100  is  designated  the  barosphere ,  or  tropopause, 
and  it  is  a  region  of  transition  above  which  very  little 
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vertical  mixing  occurs.  Above  the  barosphere,  from  100  km  to 
around  1000  km,  lies  the  heterosphere .  Due  to  the  absence  of 
mixing  in  this  region,  the  ratio  of  constituents  varies  with 
altitude.  At  extremely  high  altitudes  of  1000  kilometers  or 
so,  the  atmosphere  becomes  so  tenuous  that  particle 
trajectories  are  essentially  parabolic  orbits.  This  "outer 
atmosphere"  is  called  the  exosphere . 

The  data  analysis  in  this  paper  is  concentrated  in  the 
lower  heterosphere,  the  first  200  km  of  the  thermosphere, 
where  there  is  an  obvious  positive  temperature  gradient  and 
the  convective  mixing  is  minimal.  This  region  is  also  termed 
the  ionosphere. 

The  ionosphere  is  so  termed  due  to  the  ionization  of  the 
atmosphere  by  absorbed  solar  radiation.  The  actual  height  and 
thickness  of  the  four  layers  (D,  E,  FI ,  F2 )  vary  with  solar 
activity  and  atmospheric  composition.  The  D  region  is 
typically  from  70  km  to  90  km  during  the  day;  the  E  region  is 
typically  from  90  to  130  km;  FI  generally  exists  from  130  to 
160  km  during  the  day;  and  F2  is  above  160  km.  Each  of  these 
regions  is  formed  by  ionization  due  to  solar  radiation  at 
different  wavelengths.  For  example,  longer  wavelength 
radiation  can  penetrate  deeper  into  the  atmosphere.  Since 
solar  activity  changes  daily,  seasonally,  and  yearly,  the 
intensity  of  sunlight  at  a  particular  wavelength  and  hence  the 
depth  of  penetration  is  variable.  Ionization  depends  on  the 
composition  of  the  neutral  atmosphere  as  well.  Combining 
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these  two  factors  of  composition  and  solar  activity,  one  can 
see  that  the  flux  of  radiation  is  greatest  at  high  altitude 
but  density  is  lowest  and  vice  versa.  The  maximum  ionization 
occurs  somewhere  between  these  two  extremes,  and  the  result  is 
an  overall  layer  known  as  a  Chapman  layer.  The  composition 
and  formation  of  the  ionosphere  is  quite  complex  and  is  not 
unique  to  Earth.  An  excellent  discussion  of  planetary 
ionospheres  is  written  by  Chamberlain  (1978). 

2.  Atmospheric  Airglow 

Airglow  is  the  naturally  occurring  radiation  in  the 
Earth's  atmosphere.  The  daytime  emissions,  called  dayglow, 
are  primarily  driven  by  photoexcitation.  Also  contributing  to 
dayglow  are  emissions  which  result  from  photoelectron  impacts 
and,  to  a  lesser  degree,  photochemical  reactions.  During 
daylight  hours  photochemical  reactions  are  negligible  compared 
to  the  abundance  of  emissions  due  to  sunlight,  but  the  same 
reactions  are  significant  at  night  since  there  are  obviously 
few  direct  solar  photons  to  stimulate  transitions.  Because  the 
data  analyzed  in  this  work  were  obtained  during  the  day, 
photochemical  reactions  are  not  specifically  addressed  in  this 
work . 

The  primary  source  of  dayglow  emissions  is 
photoexcitation.  Photoexcitation  is  the  absorption  of  a 
photon  as  a  result  of  inelastic  collision  between  a  photon  and 
an  atom  or  molecule  to  raise  the  particle  to  an  excited  state. 
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When  the  photon  collides  with  the  particle,  the  photon  is 
absorbed  and  the  particle  jumps  to  a  higher  energy  level. 
When  the  particle  drops  back  to  a  lower  energy  state,  it  emits 
a  photon  with  energy  equal  to  the  transition.  This  process 
will  be  farther  described  in  the  next  section  on  atomic  and 
molecular  emissions. 

Although  not  as  important  to  dayglow  as  photoexcitation, 
photoelectron  collisions  are  significant  because  they  can 
produce  transitions  which  are  otherwise  forbidden. 
Observations  of  forbidden  transitions  can  be  used  to  find  the 
photoelectron  flux,  since  only  photoelectron  collisions  could 
excite  the  particle  to  the  upper  state.  Given  the  importance 
of  whether  a  transition  is  allowed  or  forbidden  to 
calculations  of  photoelectron  flux,  a  knowledge  of  basic 
selection  rules  is  desirable  for  the  analysis  of  dayglow 
emissions . 

It  is  also  worthy  to  note  that  some  of  the  forbidden 
transitions  which  are  prominent  in  the  low  density  unconfined 
gas  of  the  thermosphere  cannot  be  produced  in  earthbound  lab 
experiments  (Rees, 1989).  Such  emissions,  seen  in  dayglow  data 
collected  by  space-borne  equipment,  do  not  match  any  of  the 
spectra  which  have  been  observed  using  lab  or  ground-collected 
data.  These  emissions  are  identified  with  transitions  of 
known  elements,  by  wavelength  matching.  One  needs  to  be 
familiar  with  spectral  nomenclature  and  selection  rules  in 
order  to  understand  the  explanations  of  these  identifications. 
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B.  ATOMIC  AND  DIATOMIC  SPECTRA 

1.  Atomic  Emissions 

The  primary  atom-photon  interactions  are  shown  in  Figure 
2-2,  In  elastic  scattering,  no  energy  transfer  takes  place 
between  the  photon  and  the  atom  during  collision,  although  the 
photon  may  change  directions.  In  inelastic  scattering,  the 
incident  photon  may  not  only  change  direction,  but  may  impart 
energy  to  the  atom.  In  this  case,  the  atom  is  excited  to  a 
level  AE  above  its  original  state  and  the  deflected  photon 
will  have  an  energy  AE  less  than  its  incident  energy.  A 
photon  with  energy  AE  will  be  emitted  when  the  atom  relaxes  to 
the  initial  state.  Resonance  scattering  occurs  when  the 
incident  photon  is  absorbed,  raising  the  atom  to  an  excited 
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state.  A  photon  with  wavelength  equal  to  the  incident  photon 
is  emitted  when  the  atom  falls  back  to  its  initial  state. 
Fluorescence  is  the  phenomenon  which  occurs  when  the  excited 
atom  undergoes  several  downward  transitions  in  returning  to 
the  initial  state,  giving  a  cascade  of  photons.  An  atom  can 
also  return  to  the  ground  state  by  colliding  with  another 
particle,  thus  returning  to  the  ground  state  without  emitting 
a  photon.  This  is  called  collisional  de-excitation  or 
quenching . 

Similar  to  photon-atom  collisions,  electron-atom 
collisions  may  be  elastic  or  inelastic.  In  an  electron-atom 
collision,  the  scattered  electron  loses  kinetic  energy,  rather 
than  actually  being  absorbed,  when  the  atom  is  excited  to  a 
higher  energy  level.  Photon  emissions  occur  when  the  atom 
relaxes  to  the  ground  state. 

2.  Molecular  Emissions 

Molecular  spectra  may  occur  as  the  result  of  electronic 
state  changes,  but  they  may  also  result  from  purely 
vibrational  and  rotational  state  changes.  Each  of  these  types 
of  transitions  has  typical  energy  gaps,  with  electronic  being 
on  the  order  of  10  eV,  vibrational  of  .1  eV,  and  rotational  of 
.001  eV.  These  three  types  of  transitions  are  illustrated  in 
Figure  2-3.  As  seen  in  the  figure,  for  a  given  electronic 
energy  level  there  are  many  vibrational  energy  levels. 
Similarly,  for  a  given  change  in  vibrational  energy  there  are 
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Figure  2-3  Molecular  absorption  transitions: 

Transition  I  is  pure  rotational 
Transition  II,  vibrational  change  Av=l 
Transition  III,  electronic  change  of  state 
(Brehm  &  Mullin,  1989) . 

many  possible  initial  and  final  rotational  levels.  As  a 
result  of  the  many  vibrational  and  rotational  levels  which  can 
be  populated  for  any  given  electronic  transition,  the  emission 
feature  is  composed  of  many  closely  spaced  lines.  Figure  2-4 
shows  a  molecular  spectrum,  exhibiting  a  distribution  of 
emissions.  When  the  separation  of  component  lines  is  smaller 
than  instrument  resolution,  the  transition  may  appear  as  a 
hump,  as  shown  in  Figure  2.5.  This  is  the  case  for  the  data 
studied  in  this  paper. 
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The  study  of  quantum 
mechanics  explores 
transitional  possibilities, 
and  develops  a  set  of 
guidelines  for  those  emissions 
one  expects  to  observe  These 
so-called  selection  rules 
require  that  rotational 
transitions  only  occur  for 
Aj=±i,  where  j  is  the  total 
quantum  number 
( j=0 » 1 / 2 , 3 , . . . ) •  Vibrational 
transitions  occur  within  the 
confines  of  selection  rules 
Aj=±l  and  Av=±l,  where  v  is 
the  vibrational  quantum 
number  ( v=0 , 1 ,  .  .  .  j  )  .  These  rules  are  not  strictly  adhered  to, 
as  they  are  based  on  electronic  dipole  behavior,  and  some 
molecules  exhibit  quadrupole  and  magnetic  dipole 
characteristics.  There  are  also  cases  where,  though  only 
slightly  probable,  emissions  occur  when  av>±1.  Electronic 
transition  selection  rules  are  complex,  and  a  full  treatment 
is  given  to  the  development  and  explanation  of  them  by  Green 
and  Wyatt  (1965).  A  summary  of  these  is  shown  in  Table  2-1. 


Figure  2-4  Graphical 
representation  of  molecular 
spectrum . 


to  molecular  transition. 
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TABLE  2-1  SUMMARY  OF  SELECTION  RULES  FOR  MOLECULAR 
TRANSITIONS  (Green  &  Wyatt,  1965) . 


Elcctrk  Dipole  (El) 

Ma;;i-elic  Dipole  (Ml) 

Elcrtric  Quadrupole  (E2) 

Tou!  inversion 

-U  4-t  .  4-4 

+ 

Elelrnnic 

inversion 

g*-¥U 

Q  u  <-♦  U 

? 

<->  J  u*-*  u 

AJ 

0,  ±1 

0‘,±1  (0  ^  0] 

0, 

±1,  ±2) 

AN(a;,(b) 

0 

0 

0 

A2  ;a] 

(I 

0 

A. A  (a) 

o,=i 

±lif  AS  =  0, 

OifAS  =  ±1 

0.. 

±1,  ±2 

An  (a) 

±1 

0, 

±1,  ±2 

AA  (b] 

O’,  ±  1 

0,  rl 

0, 

±1,  i2 

^K  (h) 

0§,=1 

0,  =I§ 

0. 

±1§,  ±2 

‘Not  0-0.  tXot  0  -  0,  0  -  I, 

{  If  0  0  thea  i/  0.  §  Not  fo:  Z  —  I  truxiijitions. 


Derivation  of  these  rules  is  not  essential  to  this  data 
analysis,  but  an  understanding  of  the  corresponding  molecular 
nomenclature  is  necessary  in  order  to  understand  the 
discussion . 

The  component  of  total  angular  momentum  which  lies  along 
the  internuclear  axis  of  a  molecule  is  designated  A,  and  is 
used  to  label  molecular  electronic  state.  Each  state  is  given 
a  capital  Greek  letter  depending  on  its  value  of  A.  These 
designations  are  shown  in  Table  2-2.  The  spin  of  the  molecule 
is  designated  Z  (which  should  not  be  confused  with  the  E 
symbol  for  the  A=0  state) .  It  can  take  on  values  from  -S  to  S, 
where  S  is  the  resultant  of  individual  electron  spins  in  the 
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TABLE  2-2  MOLECULAR 
ELECTRONIC  STATES. 

molecule.  Multiplicity,  which  is 
2S  +  1,  is  written  as  a  left 
superscript  on  the  state  letter. 

The  ground  state  of  a  diatomic 
molecule  is  labelled  with  a 
leading  X.  Successive  excited 
states  of  the  same  multiplicity 
are  labelled  A,B,C,  and  so  forth. 

Excited  states  of  different 
multiplicity  are  alphabetically 
designated  by  small  letters.  If  these  labels  are  strictly 
followed,  then  in  general  the  selection  rules  allow 
transitions  between  like-labelled  states  (e.g. 
capital^capital) .  As  with  any  rule,  there  are  exceptions. 
For  instance,  Nj  diagrams  have  capital  and  small  letters 
assigned  conversely. 


III.  DAYGLOW  EMISSION  FEATURES 


A.  MOLECULAR  NITROGEN 

Figure  3-1  is  an  N2  energy  level  diagram  including  the 
three  transition  bands  of  interest  in  the  1800-3400A  range. 
The  states  are  indicated  by 
molecular  notation,  and  the 
transitions  are  identified  by 
name . 

1.  Vegard-Kaplan 

The  Vegard-Kaplan  system 
is  so  designated  because 
Vegard  discovered  the  system 
in  1932,  and  Kaplan  determined 
the  specific  wavelengths 
involved  a  few  years  later. 

The  emissions  range  from  125oA 
to  5325A  and  represent  the 
transition.  In  this 
transition,  the  total 
transformation  is  +-*+,  the  electronic  inversion  is  u->g,  AZ=-1, 
and  AA=0.  From  the  molecular  selection  rules  summarized  in 
Table  2-1,  this  is  a  forbidden  transition,  yet  it  has  been 
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observed  as  an  electric  dipole  transition.  A  thorough 
discussion  of  the  band  system  is  given  in  Meier  (1990). 

A-state  nitrogen  molecules  are  produced  by  cascading 
downward  from  the  B  and  C  states  or  as  a  result  of 
photoelectron  collisions.  Because  the  relationship  of 
emission  bands  to  photoelectron  flux  is  complicated  by 
cascading  N2  molecules  and  quenching  effects,  the  V-K  band  is 
not  generally  used  to  infer  photoelectron  flux. 

2.  Lyman-Birge-Hopf ield 

The  Lyman-Birge-Hopf ield  system  is  the  name  given  to  the 
a^Hg-^X^Zg*  transition  of  N^.  The  emission  bands  arising  from 
this  transition  extend  from  II60A  to  3020A  and  are  due  to  both 
electric  quadrupole  and  magnetic  dipole  transitions.  Unli>:e 
the  V-K  system,  the  upper  state  of  the  LBH  system  can  only  be 
populated  by  photoelectron  collision  excitation.  This  system 
can  therefore  be  used  to  directly  infer  the  photoelectron 
flux.  The  LBH  bands  found  in  the  data  were  particularly 
significant  in  the  wavelength  region  from  I8O0A  to  240oA.  An 
intensive  investigation  of  the  LBH  system  as  it  pertains  to 
the  MUSTANG  data  was  performed  by  Mack  (1991). 

3.  Second  Positive 

The  Nitrogen  Second  Positive  system  (C^n^-B^Ilg)  extends 
from  268oA  to  546oA.  It  is  a  particularly  important 
transition  because  the  upper  level  is  not  directly  populated 
through  photoexcitation,  but  only  by  collision.  Thus  a 
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Tneasurement  of  the  Second  Positive  intensity  profile  can  be 
used  to  infer  the  photoelectron  flux.  This  result  can  be 
compared  with  those  obtained  from  the  LBH  bands. 
Unfortunately,  there  is  no  synthetic  spectrum  model  developed 
for  the  Second  Positive.  It  was  originally  intended  to 
determine  the  N2  Second  Positive  spectrum  by  subtracting  known 
contributions  from  the  data  and  assuming  the  remaining  data 
was  Second  Positive.  This  means  of  estimating  a  spectrum  was 
found  to  be  impractical  due  to  inaccuracies  in  the  data  caused 
by  data  drops. 

B.  ATOMIC  NITROGEN 

Figure  3-2  shows  the  transitions  which  have  been  observed 
in  the  airglow  due  to  N*.  Two  transitions  are  in  the 
wavelength  range  of  this  study.  The  307oA  emission  falls  in 
the  same  range  of  wavelengths  as  the  Nj  Second  Positive 
spectrum,  so  it  could  not  be  definitively  seen  in  the  data. 
The  emission  shown  at  2143A  was  considered  in  the  synthetic 
fit  to  the  MUSTANG  data.  It  is  actually  a  doublet  at  2143. SSA 
and  2139. 68A  with  a  ratio  of  100:58,  respectively.  This 
intensity  ratio  was  determined  by  Bucsela  and  Sharp  (1989) 
from  dayglow  observation.  Theoretical  and  laboratory 
estimates  of  relative  line  strength  range  from  100:60  to 
100:23.  Since  the  data  source  for  the  Bucsela  and  Sharp  ratio 
was  similar  to  the  experimental  data  being  analyzed,  the 
100:58  relation  was  used  in  the  synthetic  model. 
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C.  ATOMIC  OXYGEN 

Figures  3-3  and  3-4  diagram  the  atomic  oxygen  emissions 
which  have  been  viewed  in  the  airglow.  Two  of  these  emissions 
fall  in  the  wavelength  range  explored  by  this  analysis. 

The  neutral  oxygen  atom  emits  a  photon  with  wavelength 
2912k  when  it  falls  from  the  ’s  state  to  the  ground  state. 
This  transition  is  depicted  in  Figure  3-3.  Singly  ionized 
oxygen,  shown  in  Figure  3-4,  produces  a  line  at  247oA  when  it 
falls  from  the  state  to  the  ‘'S  ground  state. 
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D.  NITRIC  OXIDE 


Nitric  oxide  emissions,  shown  in  Figure  3-5,  dominate  much 
of  the  dayglow  data.  There  are  four  primary  bands  in  the 
range  I8O0A  -  340oA,  7,  5,  (3,  and  e.  The  (8  bands  are  several 
orders  of  magnitude  less  intense  than  the  other  emissions,  so 
they  have  been  neglected. 

1.  Delta 

The  C^n-»X^n  transition  band  system,  designated  6,  ranges 
from  162iA  to  4288A  (Barth, 1965) .  Using  standard  notation, 
the  upper  electronic  state's  vibrational  level  is  designated 
v' ,  and  the  lower  state  vibrational  level,  v" .  The  6 
transition  band  includes  vibrational  levels  corresponding  to 
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Nitric  oxide  will  dissociate 


v'=0,l,2,3,4  and  v"=0 , 1 , . . . 23 . 
at  v'>0  for  C^n,  so  v'  is  limited  to  v'=0  and  v"<14  in 
creating  the  synthetic  6  spectrum. 

2 .  Gamma 

The  7  bands  of  nitric  oxide  extend  from  1873A  to  6126A. 
The  emissions  arise  from  the  transition,  and  the 

ground  state  has  v"=0,...23  as  previously  stated.  The 
dissociation  energy  for  nitric  oxide  is  lower  than  the  energy 
of  the  v'=4  level  in  the  A  state,  so  only  transitions  arising 
from  v'=0,l,2,3  contribute  to  the  synthetic  y  spectrum. 


Figure  3-4  Energy  level  diagram  for  O*  (Rees,  1989). 
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IV.  THE  EXPERIMENT 


A.  Introduction 

The  middle  ultraviolet  spectra  analyzed  in  this  thesis 
were  obtained  from  the  NPS  MUSTANG  spectrograph  flown  on  a 
NASA  sounding  rocket  launched  March  30,  1990,  from  White  Sands 
Missile  Range,  New  Mexico.  The  spectrograph  was  launched  on 
a  Terrier  boosted  Black  Brant,  which  carried  the  payload  to  an 
apogee  altitude  of  320  km.  The  spectrograph  was  looking  in 
the  anti-solar  direction  with  an  observation  zenith  angle  of 
90°.  From  320  to  180  km  on  the  downleg  no  useful  data  were 
obtained  due  to  the  payload  attitude. 

B.  Instrument  Description 

The  MUSTANG  (Middle  Ultraviolet  Spectrograph)  is  a 
modified  1/8  m  Ebert-Fastie  spectrograph  with  a  1/8  m 
telescope.  The  entrance  slit  is  5mm  x  140  and  a  loA 
wavelength  resolution  is  provided  by  a  1200  g/mm  grating 
ruling.  The  detector  assembly  consists  of  an  ultraviolet-to- 
visible  image  intensifier  using  a  quartz  input  window,  a  CsTe 
photocathode , and  a  fiber  optic  output  window.  The  detector  is 
an  array  of  512  photodiodes  giving  voltage  outputs 
proportional  to  intensity  using  a  50  msec  integration  period. 
This  voltage  output  is  fed  serially  to  an  A/D  converter  and 
then  to  rocket  telemetry  via  a  FIFO  shift  register. 
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Calibration/sensitivity  curves  were  obtained  using  a 
deuterium  lamp  for  the  wavelength  range  1800-2800  A  and  a  FEL 
tungsten  lamp  for  2400-3400A.  A  more  precise  calibration  of 
the  I8OO-20O0A  range,  which  is  affected  by  0^  absorption,  was 
performed  in  a  vacuum  chamber  (Mack, 1991). 

C.  Data  Collection 

Over  8000  spectra  were  collected  by  MUSTANG.  These 
spectra  were  then  averaged  into  25  altitude  bins,  disregarding 
spectra  which  were  contaminated  by  the  instrument's 
electronics.  Each  of  the  25  resultant  spectra  covers  a  10  km 
layer  and  is  identified  by  its  central  altitude. 

Because  of  a  timing  difference  between  clocking  and 
reading  functions,  only  480  of  the  512  photodiodes'  responses 
were  recorded.  It  was  first  determined  that  every  seventeenth 
word  was  dropped,  and  these  blanks  were  inserted  as  averages 
of  the  two  adjoining  data  values  (Clayton, 1990) .  This  still 
did  not  produce  a  suitable  wavelength  match  with  known 
emissions,  so  a  mathematical  spline  was  performed  on  the  data 
in  order  to  match  a  uniform  wavelength  grid  (Mack, 1991). 
Although  this  provided  an  excellent  fit  to  assumed  models,  it 
had  no  physical  justification.  Laboratory  testing  was 
performed  by  Quint  (1991)  to  determine  which  values  were  lost 
in  the  timing  error.  His  results  accounted  for  the  32 
expected  data  dropouts,  but  two  additional  dropouts  remained 
unexplained.  These  two  shifts  were  inserted  at  positions 
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which  were  consistent  with  the  findings  of  Quint  and  also 
provided  the  best  agreement  with  three  previous  studies.  It 
should  be  pointed  out  that  this  was  not  confirmed  with  the 
electronics.  The  corrected  spectra  are  delineated  in  Appendix 
A,  noting  each  dropped  pixel.  At  these  spots  an  average  of 
the  adjacent  values  was  inserted  into  the  original  480  data 
points,  to  obtain  an  array  of  512  values. 

The  adjusted  spectra  were  then  divided  by  instrument 
sensitivity.  The  first  six 
values  were  set  to  zero,  since 
the  uncertainty  in  sensitivity 
of  the  instrument  at  the 
endpoint  is  extremely  large. 

It  was  found  that  one  further 
correction  had  to  be  made  to 
the  data.  In  both  pre-launch 
and  post-launch  high  voltage 
dark  count  spectra,  two 
voltage  spikes  existed.  These 
hot  spots  occurred  at  points 
corresponding  to  2850A  and 
3373A,  as  shown  in  Figure  4-1 
These  dark  count  hot  spots 
were  converted  to  R/A  and  then  subtracted  from  the  averaged 
spectra  before  analysis. 
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Figure  4-1  MUSTANG  "hot 
pixels"  plotted  as  intensity 
(in  R/A)  versus  wavelength. 
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V. 


DATA  ANALYSIS  &  RESULTS 


A.  TECHNIQUE 

1.  Introduction 

a.  Fitting  Method 

Synthetic  spectra  were  fit  to  the  data  using  a  least 
squares  fit  process  (hereafter  referred  to  as  GRIDLS) .  It  was 
developed  by  Bevington  (1969)  and  applied  in  limited  portions 
of  the  wavelength  range  by  Clayton  (1990),  Andersen  (1990), 
and  Mack  (1991).  This  method  is  an  extension  of  a  basic 
least-squares  regression,  minimizing  error  between  the  fit  and 
the  data  as  measured  by  the  variable  with  weighting  for 
confidence  in  the  data.  The  program  fits  a  linear  combination 
of  synthetic  spectra  to  the  data,  with  each  synthetic  spectrum 
having  an  independent  coefficient,  or  scale  factor.  The 
program  is  adaptable  for  any  number  of  desired  parameters. 
The  final  analysis  used  nine  parameters,  which  will  be 
discussed  in  the  Results  section.  The  fitting  program 
(GRIDLS)  and  its  accompanying  subprograms  are  in  Appendix  B. 

The  output  of  the  program  is  the  scale  factors  for  the 
linear  combination  of  models.  What  each  models'  scale  factors 
represent  depends  on  the  normalization  dimensions  used  in 
creating  the  corresponding  synthesized  spectrum.  The  model 
spectra  do  not  all  have  the  same  dimensional  units.  In  the 
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case  of  atomic  spectra,  the  scale  factors  are  absolute 
intensity  measured  in  Rayleighs.  For  the  VK  and  LBH  bands  the 
scale  factors  are  a  measure  of  relative  intensity  and  must  be 
converted  to  absolute  values  to  obtain  an  intensity  profile. 
The  Nitric  Oxide  scale  factors  are  column  density,  in 
molecules  per  square  centimeter. 

Jb.  Synthetic  Spectra 

Previous  analysis  of  the  March  1990  MUSTANG  data  includes: 
an  analysis  of  NO  at  2000-2500A  (Clayton ,  1990)  ;  analysis  of  01 
2912k  and  Oil  2470A  emissions  (Andersen, 1990) ;  and  a  study  of 
LBH  and  NO  at  1800-2100A  (Mack, 1991) .  An  analysis  combining 
all  the  contributing  spectra  and  spanning  the  full  MUSTANG 
wavelength  range  had  not  been  performed,  but  previous  work 
provided  many  of  the  computer  programs  as  well  as  "initial 
guesses"  for  the  fitting  routine.  Other  models  and  programs 
arose  as  a  result  of  the  fitting  analysis.  Computer  models 
written  by  Cleary  (1986)  were  used  for  the  generation  of 
Nitrogen  Lyman-Birge-Hopf ield  and  Nitric  Oxide  7,  6,  and  e 
band  emissions.  The  Nitrogen  Vegard-Kaplan  model  was  written 
by  Siskind  and  Barth  (1987). 

The  synthetic  spectra  for  individual  atomic  emissions  were 
created  by  convolving  appropriate  impulse  functions  with  the 
slit  function  of  the  instrument.  As  already  discussed,  in  the 
case  of  the  N*  2143A  doublet  a  pair  of  impulse  functions  in 
the  ratio  100:58  was  used.  The  previously  identified  atomic 
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contributions  to  the  MUSTANG  spectrum  include  01  2912k,  011 
247oA,  and  Nil  2143A. 

One  other  atomic  spectrum,  a  prominent  line  emission  at 
2853A  at  higher  altitudes,  has  been  included  in  the  fitting 
routine  as  a  result  of  observation.  Elements  with  emissions 
at  that  wavelength  include  nickel,  copper,  manganese,  iron, 
sodium,  cobalt,  and  oxygen.  Of  these,  oxygen  is  the  most 
abundant  in  the  atmosphere  where  the  emission  is  seen,  above 
200km.  Doubly  ionized  oxygen  has  an  emission  at  2853A  that 
although  not  previously  identified  in  dayglow  data,  is 
classified  "predicted"  by  Kelly  (1979). 

A  collection  of  atomic  emissions  observed  during  the 
fitting  process,  but  not  verified  by  other  literature,  was 
classified  as  a  separate  spectrum.  As  discussed  in  Chapter 
II,  an  atomic  emission  is  a  single  line  vice  the  "hump"  of  a 
molecular  transition.  Even  after  convolution  with  the 
instrument  slit  function,  the  form  of  an  atomic  emission 
feature  is  distinctive  from  a  molecular  band.  There  were 
several  areas  in  the  data  which  matched  the  form  at  a 
sufficient  number  of  altitudes  as  to  suggest  an  atomic 
emission.  These  wavelengths  were  then  checked  against  tables 
for  possible  origins  and  included  in  the  synthetic  spectrum  if 
a  likely  emission  existed,  based  on  atmospheric  composition 
and  transition  probability.  The  resultant  synthetic  spectrum 
(designated  "Other  Atomics")  is  shown  in  Figure  5-1.  The 
program  used  to  create  this  spectrum,  along  with  possible 
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Figure  5-1  "Other  Atoinics"  synthetic  spectrum 
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identifications  for  each  emission,  is  in  Appendix  B.  None  of 
these  emissions  occurred  with  sufficient  intensity  and 
regularity  to  merit  being  considered  as  a  separate  parameter. 
They  are  included  only  to  facilitate  the  fitting  procedure. 

2.  General  Programming  Technique 
a.  Fitting  Routine 

The  GRIDLS  program  which  successfully  handled  72  points  of 
data  with  three  parameters  (Mack,  1991)  was  overloaded  by  the 
inclusion  of  all  512  points  of  wavelength  and  nine  parameters 
(i.e.,  the  scale  factors  for  V-K,  01,  Oil,  OIII,  LBH,  NO,  Nil, 
Other  Atomics,  and  Background) .  Several  modifications  were 
made  to  improve  the  data  handling  capabilities  of  the  program, 
and  other  characteristics  were  accepted  as  limitations. 

Previously,  the  program  would  al'ow  negative  scale 
factors,  effectively  allowing  subtraction  of  a  spectrum.  It 
also  would  enter  an  infinite  loop  if  the  scale  factor  equalled 
zero.  These  problems  were  eliminated  in  the  program  by 
requiring  that  zero  be  the  minimum  scale  factor  allowed,  and 
by  performing  a  simple  global  minimum  check  to  prevent  a  zero 
scale  factor  from  causing  an  infinite  loop. 

The  program  had  great  difficulty  converging  to  a  value  if 
the  "initial  guess"  for  the  scale  factor  was  far  from  the 
minimum  deviation  value.  This  problem  was  not  corrected,  but 
a  test-plotting  program  was  written  to  allow  one  to  visually 
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match  the  synthetic  curve  to  the  data  so  as  to  obtain 
reasonable  first  guess  values. 

The  standard  function  of  instrument  response  was  an  11- 
element  symmetric  slit  function  (Clayton,  1990;  Andersen, 
1990;  Mack,  1991)  .  When  convolved  with  a  512-element 
synthetic  spectrum,  this  slit  function  would  sufficiently 
represent  the  instrument's  response  to  the  modeled  molecular 
spectrum.  However,  in  the  case  of  an  atomic  emission,  which 
occurs  as  a  line  at  a  particular  wavelength  vie  a  band  of 
lines,  the  slit  function  could  shift  the  apparent  peak  as 
much  as  1.1 k.  The  symmetric  function  centered  the  response  on 
the  applied  pixel,  a  3.133A  range,  instead  of  at  the 
particular  desired  wavelength.  A  new  set  of  five  slit 
functions,  shown  in  Figure  5-2,  was  developed  to  correctly 
model  atomic  emissions,  allowing  for  emissions  which  are  not 
centered  on  a  pixel. 

A  major  limitation  of  the  previous  version  of  the  GRIDLS 
fitting  routine  was  that  there  existed  no  absolute  measurement 
of  "goodness  of  fit."  Although  the  program  minimized  a 
value,  the  values  were  not  weighted  and  thus  could  not  be 
compared  between  runs  or  altitudes  to  determine  goodness  of 
fit  as  parameters  were  tweaked.  Creating  a  proper  weighting 
for  the  data  evaluation  so  that  these  comparisons  could  be 
made  was  a  major  task. 
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Figure  5-2  Functions  to  simulate  instrument  slit  response. 
Delta  function  is  placed:  a)  centered  pixel  right  of 
center  c)'4  pixel  right  d)^  pixel  left  e)^  pixel  left. 
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b.  Weighting  Program 

The  GRIDLS  program  allows  for  statistical  weighting  of  the 
data,  instrumental  weighting,  or  no  weighting.  Instrumental 
weighting  was  chosen  for  its  advantage  of  stressing  datapoints 
which  have  small  uncertainties.  The  first  step  in  producing 
instrumental  data  weights  was  to  determine  the  standard 
deviations  of  the  mean  (sdom)  for  every  data  point  in  each 
altitude  bin.  A  second  source  of  instrumental  uncertainty  was 
instrument  response.  It  was  determined  that  the  detector  had 
a  linear  response  for  most  output  values,  but  reached 
saturation  above  900  decimal  units  (divided  by  instrument 
sensitivity  to  get  R/A) ,  so  data  points  in  the  saturation 
range  were  given  weights  based  on  the  sdom  being  equal  to  the 
data  value  itself.  The  exception  to  this  was  when  the  data 
value  was  so  small  as  to  give  a  large  weight  to  the  point 
(since  weight  oc  sdom’’)  ,  in  which  case  the  weight  was  assigned 
1/10  the  largest  weight  factor  for  that  altitude  (Cleary, 
private  communication,  1991)  .  Additionally,  the  data  points 
which  were  merely  averages  of  adjoining  values  due  to  lost 
data  points  were  given  weights  based  on  the  sdom  being  equal 
to  the  data  value.  Since  the  sdom  and  saturation  points  are 
different  for  each  altitude,  a  set  of  values  for  weighting  of 
the  512  data  points  was  stored  for  each  data  spectrum. 

The  data  pixels  corresponding  to  nitrogen  Second  Positive 
emissions  had  no  synthetic  input  being  matched  to  them.  The 
pixel  corresponding  to  each  wavelength  of  these  emissions  as 
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determined  by  Barth  (1965),  plus  a  buffer  of  ±1  pixel  to 
account  for  instrument  response,  were  given  essentially  zero 
weight.  The  final  weight  program,  designated  WEIGHTS2,  is 
Appendix  C. 

B.  RESULTS 

1.  Comparison  of  Data  and  Fit 

a.  Nitric  Oxide  Bands 

In  the  plots  of  the  fits,  nitric  oxide  is  viewed  as  one 
entity  rather  than  plotting  the  y,  S,  and  e  bands  separately. 
This  is  for  ease  in  reading  the  plots.  Figure  5-3  shows  the 
components  of  the  nitric  oxide  synthetic  spectrum  for  115km. 
The  bands  are  not  equal  contributors,  as  the  6  band  has  a  25% 
fluorescent  efficiency  factor. 

Figure  5-4  is  a  plot  of  the  observed  data  and  its 
corresponding  fit.  The  data  is  indicated  by  the  solid  line, 
and  the  fit  obtained  from  the  computer  procedure  by  a  dotted 
line.  The  dashed  line  is  that  portion  of  the  synthetic  fit 
which  is  nitric  oxide.  From  the  plots  in  Figure  5-4,  one  can 
see  the  importance  of  nitric  oxide  with  decreasing  altitude. 
In  the  105km  plot,  nitric  oxide  accounts  for  almost  all  of  the 
dayglow  from  1800-2600A,  but  at  295km  the  NO  is  less 
prevalent . 

b.  Nitrogen  Bands 

Figure  5-5  shows  an  intensity  plot  (in  R/A)  of  the  LBH 
contribution  to  the  fit,  along  with  the  observed  spectrum,  for 
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Figure  5-5  Cor,parison  of  data  &  LBH  synthetic  spectrum  for 
295,  155,  105km. 
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three  different  altitudes.  From  these  plots,  one  can  see  that 
the  Lyman-Birge-Hopf ield  band  of  N2  emissions  dominates  the 
wavelength  range  1800-2000A. 

Figure  5-6  shows  a  similar  altitude  progression  for  the 
Vegard-Kaplan  system,  emphasizing  the  (0,5)  transition  at 
26O5A,  VK  emissions  are  prominent  from  2200-3300A,  and  in 
general  all  the  emission  features  show  the  same  characteristic 
growth  as  exhibited  by  the  (0,5)  transition, 
c.  Oxygen  Atomic  Emissions 

The  01  2972A  emission  as  fit  to  the  data  at  three 

different  altitudes  is  shown  in  Figure  5-7.  The  poor  fit  of 
this  feature  is  addressed  in  the  Discussion  section.  The  Oil 
247oA  1  ine,  shown  in  Figure  5-8,  becomes  more  important  with 
increasing  altitude. 

2.  Profiles 

a.  Nitric  Oxide 

The  column  density  profile  for  Nitric  Oxide  is  shown  in 
Figure  5-9,  with  the  corresponding  values  given  in  Table  5-1. 
Values  obtained  from  upleg  data  are  indicated  with  squares  and 
downleg  values  by  triangles.  There  is  a  systematic 
uncertainty  of  15%.  The  values  obtained  show  general 
agreement  in  both  magnitude  and  slope  with  the  column 
densities  obtained  hy  Mar'W  (19P1)  ,  Clayton  (1990)  and  Andersen 
(1990),  all  of  whom  analyzed  subsets  of  the  MUSTANG  data.  The 
maximum  and  minimum  values  best  agreed  with  those  obtained  by 
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Column  Density  Profile  of  Nitric  Oxide 
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Figure  b-9  Nitric  oxide  coluirn  density  profile. 
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Cofumn  Density  (mo(ecules/crTr2) 


TABLE  5-1  NO  COLUMN  DENSITIES 

Mack  (1991),  with  the 
intermediate  values  best 
matching  the  fit  obtained  by 
Cla/ton  ( 1990) .  The 

significance  of  this  agreement 
is  that  a  full  spectrum 
fitting  procedure  (i.e., 
including  all  species  and  all 
wavelengths)  is  not  required 
to  accurately  determine  the 
Nitric  Oxide  density. 

b.  Nitrogen  Bands 
The  scale  factors  for  LBH 
bands  were  converted  to 
absolute  intensity  using  the 
algorithm  and  computer 
programs  in  Appendix  D.  The  resultant  LBH  intensity  profile 
is  shown  in  Figure  5-10,  and  is  a  measuremenc  of  total  LBH 
intensity  as  it  includes  all  a-*X  transitions.  The  intensity 
profile  for  VK  can  only  be  calculated  for  a  particular  v'—v" 
band,  and  Figure  5-11  shows  this  calculation  for  one  of  the 
prominent  systems,  the  (0,5). 
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Figure  5-10  Nitrogen  LBH  intensity  profile. 
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Intensity  (Rayleighs) 
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Figure  5-12  Intensity  profiles  for  01,  Oil,  OIII,  Nil 
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c.  Atomic  Emissions 


The  intensity  profiles  for  01  2972A,  Oil  247oA,  OIII 

2853A,  and  Nil  2143A 

are  shown  in  Figure  TABLE  5-2  ATOMIC  INTENSITY  VALUES 
5-12,  with  the  values 
given  in  Table  5-2. 

The  01  2972A  profile 
shows  a  distinctive 
peak  at  about  150kin; 
the  O'"  247oA  profile 
shows  a  broad  peak  at 
about  250km;  the  OIII 
profile  suggests  a 
peak  at  around  300km; 
and  the  Nil  reaches 
its  maximum  around 
150km. 

The  01  2912k  line 
emission  profile  was 
found  to  have  very  little  agreement  with  previous  analysis. 
The  maximum  found  by  this  analysis  was  50%  of  that  found  when 
fitting  only  the  OI  2912k  line  (Andersen,  1930).  This  is  due 
to  the  uncertainty  of  the  wavelength  for  the  data,  since  two 
of  the  data  dropouts  were  in  this  region.  The  overall 
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uncertainty  in  the  01  profile  should  reflect  the  data  drops, 
thus  caution  is  urged  for  any  subsequent  analysis  of  this 
emission  profile. 

3.  Discussion 

a.  V-K  Quenching 

Figure  5-13  shows  an  overlay  of  the  nitrogen  LBH  and  VK 
profiles,  with  the  VK  scaled  up  by  a  factor  of  9.5  for 
comparison.  The  most  interesting  aspect  of  the  nitrogen 
profiles  is  that  they  have  roughly  the  same  slope  above  200 
kilometers,  but  the  VK  intensity  drops  off  rapidly  below  this 
height.  This  is  due  to  the  quenching  effect  of  N2,  NO,  O,  O^, 
and  N  at  the  lower  altitudes.  From  auroral  studies,  atomic 
oxygen  is  thought  to  be  the  primary  quenching  ^gent  with 
molecular  oxygen  playing  a  role  below  135km  (Sharp,  1971)  .  As 
previously  discussed,  the  upper  state  for  the  VK  emissions  is 
the  Nj  A-state.  The  A-state  to  ground-state  energy  for  N2 
corresponds  to  an  excitation  energy  for  atomic  oxygen.  Below 
200  km,  the  O  density  is  high  enough  that  during  their 
metastable  lifetime  in  the  excited  state,  N2  molecules  collide 
with  and  excite  oxygen,  thus  losing  their  energy  without 
emitting  a  photon.  This  is  believed  to  be  the  first 
simultaneous  observation  of  these  two  profiles  in  the  dayglow. 

b.  Adjustments  to  Franck-Condon  Factors 

The  first  set  of  spectra  to  be  fitted  to  the  data  included 
Nitric  Oxide  7,  5,  and  e  bands.  Nitrogen  LBH  and  V-K  bands,  01 
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Intensity  Profiles  for  LBH  &  VK 
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Figure  5-13  Intensity  plots  for  LBH  &  VK  showing  quenching 
of  VK, 
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Intensity  (Rayleighs) 


2912k,  011  247oA,  and  OIII  2853A  emissions.  The  first  fine- 
tuning  of  the  fit  was  accomplished  by  plotting  the  fit  and 
data  simultaneously  in  40oA  segments  and  looking  at  each  band 
for  trends  in  undershooting  or  overshooting  at  a  particular 
altitude,  or  for  all  altitudes  at  a  particular  wavelength. 
Any  trend  was  then  adjusted  up  or  down  in  the  original  model 
spectrum,  a  procedure  equivalent  to  varying  the  Franck-Condon 
factor  (the  factor  which  quantifies  probability  of  transition 
for  the  synthetic  models) .  This  process  was  very  tedious  and 
time  intensive,  but  a  second  fitting  yielded  an  average 
decrease  of  30%  in  x^- 

After  the  first  set  of  changes  had  been  made  to  the 
Franck-Condon  factors,  the  same  process  of  plotting  and 
observing  trends  was  repeated.  Several  unidentified  peaks 
were  seen  in  the  data  which  are  believed  to  be  unidentified 
atomic  emissions,  since  their  profiles  matched  that  of  an 
impulse  function  convolved  with  the  instrument  slit  function. 
Other  than  the  well-documented  Nil  2143A  doublet,  which  was 
added  as  a  synthetic  spectrum,  these  emissions  had  not  been 
previously  reported  in  the  literature,  so  the  "Other  Atomics" 
spectrum  was  created  (see  discussion  above) .  The  two 
parameters.  Nil  and  "Other",  were  included  in  the  fitting 
routine.  Once  these  parameters  were  included,  the  adjustments 
to  the  Franck-Condon  factors  found  for  the  Nitric  Oxide  bands' 
synthetic  spectra  confirmed  those  used  by  Mack  (1991)  and 
Clayton  (1990).  The  adjustments  to  the  Franck-Condon  factors 
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TABLE  5-3 


are  summarized  in  Table  5-3.  These  changes  to  the  accepted 
probabilities  of  transition  are  believed  to  be  actual 
corrections  and  not  just  a  function  of  the  instrumentation. 
It  is  expected  that  further  independent  measurements  will 
verify  these  corrections. 

The  intensity  profiles  were  derived  from  the  final 
parameters  of  the  fitting  routine,  but  analysis  of  the  data 
includes  looking  not  only  at  the  profiles  but  at  the  fits 
themselves.  A  complete  set  of  plots  of  the  data  showing  a 
comparison  of  data  and  fits,  as  well  as  accompanying  plots 
showing  the  breakdown  of  the  fits  by  constituent,  are  in 
Appendix  E.  Figures  5-14  and  5-15  show  sample  fits  of  the 
entire  wavelength  range  at  low  altitude  and  high  altitude.  It 
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is  apparent  from  these  figures  that  the  overall  fit  was 
excellent . 

It  should  be  pointed  out  that  there  is  a  disagreement 
between  the  data  and  fit  at  the  wavelengths  where  the  nitrogen 
Second  Positive  emissions  are  present-  Additionally,  the  only 
other  poorly  fit  feature  was  the  293oA  to  299oA  range  at 
altitudes  above  200  km.  This  is  one  of  the  areas  one  would 
like  to  be  able  to  fit  well,  since  it  includes  the  01  2912k 
emission.  Unfortunately,  two  of  the  dropped  data  points  are 
at  2969. ik  and  2912.9k,  exactly  where  the  peak  of  01  is 
expected.  Also,  there  seems  to  be  another  emission  at  about 
2955A  which  is  not  explained  by  the  VK  bands  or  01  emission, 
creating  three  peaks  in  the  data  where  theory  predicts  only 
two.  This  feature  is  shown  in  Figure  5-16. 


Wovelength(A) 


Figure  5-16  Averaged  spectrum  for  285km.  The  peak  at  2940A 
is  VK(0,7),  the_ peak  at  2980A  is  due  to  01  2912k,  and  the 
feature  at  2955A  is  unidentified. 
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c.  Temperature  Profile 

In  addition  to  the  column  density  profile,  intensity 
profiles,  and  fits  of  model  spectra  to  data,  a  temperature 
profile  also  was  determined.  The  temperatures  described  by 
Mack  (1991)  were  used  as  initial  guesses.  Several  of  the 
temperatures  at  the  lowest  altitudes  were  increased  during 
fitting  to  improve  the  match  of  synthetic  spectra  to  data. 
The  resultant  temperature  profile  is  shown  in  Figure  5-17. 
The  temperatures  determined  by  the  fitting  procedure  are  all 
considerably  hotter  than  predicted  by  MSIS.  In  particular, 
the  asymptotic  temperature  is  400°  hotter.  This  is  possibly 
due  to  a  magnetic  storm  which  occurred  about  three  hours  prior 
to  launch. 


Tempercture  (Kelvin) 


Figure  5-17  The  temperatures  used  by  the  analysis  are 
indicated  by  symbols.  The  solid  line  is  the  MSIS 
temperature  profile  for  the  day  of  the  launch. 
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VI .  CONCLUSIONS 


A.  SUMMARY  OF  FINDINGS 

Thermospheric  dayglow  data  from  a  rocket-borne 


spectrograph 

were 

analyzed  in 

this  study. 

Computer 

synthesized 

spectra 

were  created 

for  the  Nj 

Lyman-Birge- 

Hopfield  and 

Vegard- 

•Kaplan  bands; 

NO  7,  S,  and 

€  bands;  01 

2912k,  on  247oA,  and  OIII  2853A  atomic  lines;  and  Nil  2143A 
doublet  feature.  A  spectrum  of  other  possible  atomic 
emissions  and  a  background  factor  were  also  considered  in  the 
fit.  These  synthetic  models  were  fitted  to  the  observed  data, 
then  the  scaling  factors  were  used  to  determine  an  intensity 
profile  (column  density  profile  in  the  case  of  nitric  oxide) . 

The  intensity  and  column  density  profiles  obtained  were 
within  10-20%  of  the  values  acquired  by  simpler  fitting 
routines  in  all  cases  except  the  01  2912k  emission.  Thus  the 
analysis  supports  the  concept  of  fitting  of  just  one  feature 
to  arrive  at  density  and  intensity  profiles.  A  relationship 
was  noted  between  VK  and  LBH  emissions  that  would  seem  to 
support  the  theory  of  quenching  of  the  nitrogen  VK  system. 
The  data  above  roughly  300oA  could  not  be  modelled  without  a 
Nj  Second  Positive  synthetic  spectrum.  The  loss  of  datapoints 
was  a  major  detriment  in  fitting  the  01  2912k  emission,  as 
there  were  two  data  dropouts  where  a  peak  was  expected.  The 
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temperature  profile  which  best  fit  the  data  was  higher  at  low 
altitudes  than  previous  analysis  had  predicted. 

B.  TOPICS  FOR  FURTHER  INVESTIGATION 

With  another  launch  scheduled  for  March  1992,  further 
spectral  analysis  of  the  March  1990  data  is  imprudent. 
Contemplating  the  prospective  data,  a  synthetic  spectrum  for 
the  N2  Second  Positive  would  enhance  the  analysis.  It  is 
expected  that  there  will  be  no  problem  with  data  dropouts  on 
the  next  flight  (Quint,  1991) ,  so  a  repeat  of  this  analysis 
(with  the  new  data)  will  greatly  improve  the  accuracy  of 
fitting  of  the  01  emission  and  possibly  show  that  it  can  be 
accurately  fitted  over  a  limited  wavelength  range. 

The  quenching  of  VK  below  200km  should  be  paralleled  by 
the  density  of  the  quenching  agent,  as  well  as  a  rise  in 
emissions  corresponding  to  the  relaxation  of  the  quenching 
agent  from  its  excited  state.  The  quenching  phenomenon  showed 
promise  for  rewarding  study,  but  was  not  the  subject  of  this 
thesis . 

Looking  toward  the  future,  the  MUSTANG  will  be  part  of  an 
operational  satellite  in  the  late  1990's,  providing  a  constant 
flow  of  data  for  analysis.  Considering  the  numerous  new 
mathematical  theories  for  computer  network  calculations,  there 
may  be  a  method  better  suited  for  fitting  large  quantities  of 
data  than  the  GRIDLS  routine.  A  modification  to  the  fitting 
programs  to  place  less  stringent  requirements  on  the  initial 
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programs  to  place  less  stringent  requirements  on  the  initial 
guess  and  to  provide  a  more  formal  check  for  global  minimum 
would  greatly  improve  the  effectiveness  of  the  existing 
routine.  Additionally,  automation  of  the  analysis,  including 
derivation  of  the  electron  density  profile,  would  certainly 
add  to  the  viability  of  MUSTANG  as  a  provider  of  geophysical 
parameters  for  national  defense  agencies. 
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APPENDIX  A 


The  following  table  shows  the  unshifted  and  shifted  array 
elements  of  data,  and  their  corresponding  wavelengths. 


Record:  [mustang] iii90  avg/um90  avg 

Entry  M 
(pixel) 
Shitted 

Entry  # 
(pixel) 
Unshifted 

Wavelength 

A 

0 

0 

1,798.00 

1 

1 

1,801.13 

2 

2 

1,804.27 

3 

3 

1,807.40 

4 

4 

1,810.53 

5 

5 

1,813.67 

6 

6 

1,816.80 

7 

7 

1 ,819.93 

8 

8 

1,823.06 

9 

9 

1,826.20 

10 

10 

1,829.33 

11 

11 

1,832.46 

12 

12 

1,835.60 

13 

13 

1,838.73 

14 

14 

1,841.86 

15 

15 

1,845.00 

16 

1,848.13 

17 

16 

1,851.26 

18 

17 

1,854.39 

19 

18 

1,857.53 

20 

19 

1,860.66 

21 

20 

1,863.79 

22 

21 

1,866.93 

23 

22 

1,870.06 

24 

23 

1,873.19 

25 

24 

1,876.33 

26 

25 

1,879.46 

27 

26 

1,882.59 

58 


59 


60 


61 


63 


64 


65 


66 


67 


68 


59 


70 


71 


503 


470 


504 

471 

3,377.032 

505 

472 

3,380.165 

506 

1,71 

3,383.298 

507 

474 

3,386.431 

508 

475 

3,389.564 

509 

476 

3.392.697 

510 

477 

3,395.830 

511 

3,398.963 
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APPENDIX  B 


Appendix  B  contains  the  least  squares  fitting  program  and 
its  subroutines.  Additionally,  the  synthetic  spectrum 
generation  routines  are  given. 
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; Walden , 1991 ,  modified  from  Mack, 1991 

@getfunc 

@f chi sqr 

pro  nimpfit_all 

;  This  program  fits  data  for  the  1800  to  340 OA  range  with  models  of  the 

;  NO  gamma,  delta,  and  epsilon  bands,  N2  LBH  emissions,  N2  Vegard- 

;  Kaplan  bands,  oxygen  emissions  at  2470,  2972,  &  2853a,  N+  doublet 

;  at  2143A,  various  atomic  emissions,  and  the  background  level 

;  through  an  iterative  process.  It  uses  gridls.pro,  getfunc.pro,  and 

;  fchisqr.pro  as  subprograms;  it  uses  parameter.dat  as  an  in/out  file; 

;  and  it  writes  the  fit  by  altitude  to  nimpfit_all.dat.  The  subprograms 

;  utilize  synthetic  spectra  which  are  stored  in  allvk.dat,  all  Ibh.dat, 

;  nii_214  3doub.dat,  nodlt.dat,  nogamm.ack.dat,  noeps.dat,  oiii_58  53.dat 

;  oii_2470.dat,  oi_2972.dat,  and  otherspec.dat. 


parameter.dat  -  file  of  parameters-  records  have  25  elements  each 
record#-  contents 

0  -  temperature 

1  -  altitude  record  number  (1  through  25;1-17  up, 18-75  down) 

2  -  chi  squared  of  the  fit 

3  -  scale  factor  for  0  2972A  emission,  in  Rayle ighs ( intensity ) 

4  -  sigma  for  0  scale  factor 

5  -  scale  factor  for  0+  2470A  emission,  in  Rayleighs 

6  -  sigma  for  0+  scale  factor 

7  -  scale  factor  for  0++  2853A  emission,  in  Rayleighs 

8  -  sigma  for  0++  scale  factor 

9  -  scale  factor  for  the  N2  LBH  bands 

10  -  sigma  for  LBH  scale  factor 

11  -  scale  factor  for  the  NO  gamma,  delta,  and  epsilon  bands,  in 

molecules  per  square  centimeter  (column  density) 

12  -  sigma  for  NO  scale  factor 

13  -  scale  factor  for  the  N2  V-K  bands 

14  -  sigma  for  VK  scale  factor 

15  -  scale  factor  for  N+  2143A  emission,  in  Rayleighs 

16  -  sigma  for  N+  doublet  scale  factor 

17  -  scale  factor  for  synthetic  atomics  spectrum,  in  Rayleighs 

18  -  sigma  for  atomic  spectrum 

19  -  scale  factor  for  the  background,  in  Rayleighs  per  Angstrom 

20  -  sigma  for  background  scale  factor 

21  -  altitudes  in  kilometers 

para  -  a r r ay ( ( 2 *nte rms ) +2 )  in  same  order  as  record  order  in  file 

a  -  array ( nterms+1 )  of  temp  followed  by  scale  factors 

yfit  -  512  element  array  of  best  fit 

start-  first  record  number 

stop  -  last  record  number 

open r , 3 0 , ' ( mus tang ] m9 0_i nt . da t '  ;open  data  file 

aa  =  assoc( 30, fltarr( 512  )  ) 

openu , 3 1 , ' nimpf i t_al 1 . da t '  ;open  file  to  save  yfit 

bb=assoc(31,fltarr(512)  ) 

openu , 32 , ' pa ramete r . dat '  ;open  parameter  file 

pp=assoc( 32,fltarr(25)  ) 

openr,33, ' [walden.datafix] m90_domwt . dat ' 
dd  =  assoc( 33,fltarr(512)  ) 
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star t=0 
stop=511 

npts=( stop-start )+l 
nterms  =  9 

para=fltarr(2*nterms+3) 
a=fltatr(nterms+l ) 


;#  of  points  to  be  fit 

;#  of  terms  to  consider,  change  if 

;pararaeters  are  added 


;fill  para  with  values  for  alt 


;assign  first  element  in  a  temp 
;fill  a  with  scale  factors 


for  i=0,24  do  begin 

for  j =0 , ( 2*nte rms ) +2  do  begin 
hold  =  pp ( j  ) 
para ( j ) =hold ( i  ) 

endf or 
alt=para ( 1 ) 
count=0 
a ( 0 ) =pa  ra ( 0  ) 

for  j = 3 , ( 2 *nte rms ) +1 , +2  do  begin 
count=count+l 
a ( count ) =pa r a ( j  ) 
endf or 

y=aa( start:stop,alt) 

f 

;deltaa  is  step  size  for  estimations,  cannot  be  zero,  first  element  is  dummy 
h=a 

for  j=l, nterms  do  begin 
h( j )=a ( j  ) 

if  (h(j)  eq  0.0)  then  h(j)=1000.00 
endf or 

del taa=h/l 0 . 


mode=l 

we i ghts=dd (alt) 


;use  sdom  for  weighting 


gridls,y, we ights,npts, nterms, mode , a ,del taa , sigma a ,yfit,chisqr, start, stop 
bb( alt )=yf it 

count=0  ; save  parameters  from  best  fit 

pa  ra ( 2 ) =chi sqr 

for  j=3 , ( 2*nterms ) +1 , +2  do  begin 
count=count+l 
para ( j  )  =a { count ) 
para( j+1 )=sigmaa( count) 
endf o  r 

for  j=0 , ( 2*nte rms ) +2  do  begin 
hold=pp( j ) 
hold ( i  )  =pa ra ( j  ) 
pp( j )=hold 
endf or 


;put  parameters  in  appropriate 
;record  in  parameter.dat 


print,  'alt=’ ,alt, 'a=' ,a, 'chisqr=' ,chisqr 
endf o  r 


close , 30 
close  ,  31 
close , 32 
close , 3  3 


etu  rn 
?nd 

^gr idl s 
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;Walden,1991  modified  from  Mack, 1991 
function  getfunc,a 

;  This  function  computes  the  synthetic  spectrum  from  1800  to  3400  A 

;  for  given  scaling  factors. 

;  INPUTS 

;  a-9-Darrayof  parameters 

;  a(0)  -  temperature  at  altitude  of  interest 

;  a(l)  -  scale  factor  for  0  emission  @  2972 

;  a(2)  -  scale  factor  for  0+  emission  @  2470 

;  a(3)  -  scale  factor  for  0++  emission  @  2853 

;  a(4)  -  scale  factor  for  N2  LBH  bands 

;  a(5)  -  scale  factor  for  the  sum  of  NO  gamma,  delta, 

;  and  epsilon  bands 

;  a(6)  -  scale  factor  for  n2  Vegard-Kaplan  bands 

;  a(7)  -  scale  factor  for  N2  doublet  emission  @2143 

;  a ( 8 )  -  scale  factor  for  "Other  Atomics"  spectrum 

;  a(9)  -  scale  factor  for  background 

;  GETfunc  returns  an  array  of  511  points  to  be  compared  to  the 

;  experimentally  recorded  profile  from  1800  to  3400  A. 

;  LBH  -  The  N2  LBH  bands  convolved  with  the  slit  function 

;  NOSYN  -  The  sum  of  the  NO  gamma,  delta,  and  epsilon  bands  convolved 

;  with  the  slit  function 

;  VK  -  The  N2  V-K  bands  convolved  with  the  slit  function 

;  01  -  The  01  emission  @  2972  A  convolved  with  the  slit  function 

;  Oil  -  The  Oil  emission  @  2470  A  convolved  with  the  slit  function 

;  OIII  -  The  OIII  emission  @  2853  A  convolved  with  the  slit  function 

;  n2  -  The  N2  emission  @2143  A  convolved  with  the  slit  function 

:  AT  -  The  spectrum  of  possible  atomic  emissions 

openr , 20 , ' ALL  LBH. DAT' 
rr=assoc(20,fTtatr(512) ) 

LBH*rr(a( 0)/100) 

openr ,21, 'NOGAMMACK.DAT' 
ss=assoc(21,fltarr(512) ) 

NOGAM=ss(a( 0 )/100 ) 

openr , 22 , 'NODLT.DAT' 
tt=assoc(22,fltarr(512) ) 

NODLT=tt ( a( 0 )/100 ) 

openr ,23, 'NOEPS.DAT' 
uu=assoc(23,fltarr(512) ) 

NOEPS=uu(a( 0)/100) 

openr , 24 , 'ALLVK.DAT' 
vv=assoc(24,fltarr(512) ) 

VK  =  vv ( a ( 0 ) /1 00 ) 

openr , 25, '01  2972.DAT' 
wv.  =  assoc(  25,7ltarr(  512  )  ) 

OI=ww( 0 ) 

openr , 26 , 'Oil  2470.DAT' 
xx=assoc(26,fltarr(512) ) 

Oi  I  =  xx ( 0 ) 

openr, 27, 'OIII_2853.DAT' 
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yy=assoc(27,fltarr(512)  ) 

Olll=yy( 0 ) 

openr , 28 , 'Nil  2143DOUB.DAT' 
zz  =  assoc(28,fTtarr(512)  ) 

N2  =  zz ( 1  ) 

openr, 29, 'otherspec.dat' 
atsp=assoc(29,fltarr(512)  ) 

AT=atsp ( 1 ) 

NOSYN  =  NOGAM  +  .25*NODLT  +  NOEPS  ;25%  flourescent  efficiency  of 

;delta  band 

getfunc  =  a(4)*LBH  +  a ( 5 ) *NOSYN/3 . 1 3 3e6  +  a ( 6 ) *VK  +  a(l)*OI  $ 

+  a(2)*0II  +  a(3)*OIII  +  a(7)*N2  +a(8)*AT  +  a(9)/3.133 
;background  is  divided  by  wavelength  bin  so  that 
;scale  factor  is  absolute  intensity,  NOSYN  is  divided 
;by  3.133e6  for  same  reason 

close , 20 
close , 21 
close , 22 
close  ,  23 
close , 24 
close , 25 
close , 26 
close , 27 
close , 28 
close , 29 

return , getfunc 
end 
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; Mack , 1991 

: unction  fchisqr, y, we ights,npts, free, mode, yfit, start, stop 
This  function  evaluates  the  reduced  chi  squared 
;  for  a  fit  to  the  data. 

;  Adapted  from  Bevington,  pg  194. 

;  INPUTS 

;  y  -  array  of  data  points 

;  weights-  array  of  weightings  based  on  sdom's  of  data  points 

;  npts  -  number  of  data  points  being  fit 

;  free  -  number  of  degrees  of  freedom 

;  mode  -  determines  method  of  weighting  the  fit 

;  +1  instrumental  weighting( sdom  weights) 

;  0  no  weighting 

;  -1  statistical  weighting 

;  yfit  -  array  of  fitted  points 

;  start  -  first  record  element  to  be  fit 

;  stop  -  last  record  element  to  be  fit 

count=0 
chi sqr=0 . 

if  mode  eq  0  then  wt=l. 
if  mode  gt  0  then  wt=weights 
if  mode  It  0  then  wt=l./abs(y) 

for  i=start,stop  do  begin 

chi sqr=chi sqr  +  wt { i )*( y ( count ) -yfi t ! count ))' 2 
count=count  +  1 

endf or 

fchisqr=chisqr/free 
return, fchisqr 

?nd 
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, -Walden ,1991 

pro  Ibhcal c , wl , int , si  i  t 

openu,10, 'all_ibh.dat' 
lbh=assoc(10,fltarr(512)  ) 


Calculates  synthetic  Ibh  spectra  at  temps 
from  100  to  2000  deg  for  wavelengths  from 
1800  to  3400  A.  Results  are  convolved  with 
MUSTANG  instrument  slit  function. 


wl  =  i ndgen ( 51 2  )  *  3 . 1 3 3  +  1798 


openr,ll, ' [waldenlslit.dat' 
sl  =  assoc(li,fltarr(ll)  ) 
si i t  =  sl ( 1 ) 


i  =  0 

for  t=l 0 0 , 2 00 0 , 1 0 0  do  begin 
i  =  i  +  l 
int=wl *  0 
intl=wl *  0 
int2=wl *  0 
i nt  3=wl *  0 
int4=wl *0 
int5=wl *0 
int6=wl *  0 
i n 1 7 =wl *  0 
i nt8=wl *  0 
i nt9=wl *  0 
intl0=wl*0 
lbhsyn,t,4,9,wl,int 
lbhsyn,t,0,6,wl,int 
lbhsyn,t,5,10,wi, int 
lbhsyn,t,l,7,wl,intl0 
int=int+intl0*l .7 
Ibhsyn ,t,6,ll,wL,int 
lbhsyn,t,2,8,wl,int 
Ibhsyn, t, 3,9,wl , inti 
int=int+intl*1.203 
lbhsyn,t,4,10,wl,int 
Ibhsyn, t,0, 7, wl, int 
lbhsyn,t,5,ll,wl,int2 
i nt= int+ int2*  1  .  4 
Ibhsyn, t,l, 8, wl,int3 
int=int  +  int3*l  .  3 
Ibhsyn, t, 6, 12, wl, int 
Ibhsyn, t, 2, 9, wl, int 
lbhsyn,t,3,10,wi,int4 
int=int+int4* . 8 
lbhsyn,t,4,ll,wl,int5 
int=int+int5* . 9 
Ibhsyn, t,0, 8, wl, int 
Ibhsyn, t, 5,12,wl , int6 
int=int+int6*.757 
lbhsyn,t,l,9,wl,int 
Ibhsyn, t, 6, 13, wl, int 
Ibhsyn, t, 2, 10, wl, int 
lbhsyn,t,3,ll,wl,int7 
int  =  int  +  int7  * . 8 
Ibhsyn , t , 4 , 1 2  wl,int8 
int  =  int  +  int8*  .  9 
lbhsyn,t,0,9,wl,int 
Ibhsyn , t , 5 , 1 3 , wl , i nt  9 


Wavelength 

variables  for  tweaking  individual  transitions 


1801 . 0 
1804.6 
1817.9 
1820.8 

1835.0 

1837.2 

1853.8 


870 . 8 
877 . 7 

1887.9 

1894 . 2 

1905.4 

1910.9 

1927.8 

1945.0 

1955.9 

1962.4 

1972.6 

1980.1 

1989.6 
2006.8 

2024 . 2 

2040 
2  0  41.9 
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int  =  int  +  int9*l .  2 
Ibhsyn, t , 1 , 10 , wl , int 
lbhsyn,t,6,14,wl,int 
Ibhsyn, t, 2, 11, wl, int 
Ibhsyn, t , 3 , 12 , wl , int 
Ibhsyn, t, 4, 13, wl, int 
Ibhsyn , t , 5 , 1 4 , wl , int 

Ibhsyn, t,0, 10, wl, int 
Ibhsyn, t, 6, 15, wl, int 
Ibhsyn, t ,1,11 ,wl , int 
Ibhsyn, t, 2, 12, wl, int 
Ibhsyn, t, 3 , 13 ,wl , int 
Ibhsyn, t, 4, 14, wl, int 
Ibhsyn, t , 5 , 15 ,wl , int 
Ibhsyn, t,0, 11, wl, int 
Ibhsyn, t, 6, 16, wl, int 
Ibhsyn, t,l, 12, wl, int 
Ibhsyn, t , 2 , 13 , wl , int 
Ibhsyn, t, 3, 14, wl, int 
Ibhsyn, t, 4, 15, wl, int 
Ibhsyn, t, 5, 16, wl, int 
Ibhsyn, t, 6, 17, wl, int 
Ibhsyn, t,l, 13, wl, int 
Ibhsyn, t, 2, 14, wl, int 
Ibhsyn, t, 3, 15, wl, int 
Ibhsyn, t , 4 , 16 , wl , int 
Ibhsyn, t, 5, 17, wl, int 
Ibhsyn, t , 6 , 18 ,wl , int 
Ibhsyn, t, 2, 15, wl, int 
Ibhsyn, t , 3 , 16 , wl , int 
Ibhsyn, t , 4 , 17 ,wl , int 
Ibhsyn, t , 5 , 18 , wl , int 
Ibhsyn, t , 6 , 19 , wl , int 
Ibhsyn, t , 4 , 18 ,wi , int 
Ibhsyn, t , 5 , 19 ,wl , int 
Ibhsyn, t , 6 , 20 , wl , int 
Ibhsyn, t, 6, 21, wl, int 

print , total ( int ) 
spec=convol(int,slit) 
lbh{i)=spec/total{spec) 
endf  or 

close , 10 
close  ,  1 1 

eturn 

nd 

[ SYNSPEC ] LBHSYN 


2056.6 

2059.8 

2073.8 
2091.2 

2108.8 

2126.7 
2130 
2144 . 0 
2147 
2164 
2181 . 1 

2198.7 
2217 
2227 
2235 
2244 
2261 
2279 
2297 
2314 
2332 
2349 
2366 
2384 
2401 
2419 
2437 
2479 
2496 
2514 
2531 
2549 
2636 
2653 
2670 
2801 
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; Walden ,1991 

pro  allvkcalc ,wl , int , slit 
openr,l, ' [WALDENjslit.dat' 
openu, 3 , ' allvk.dat' 
sl=assoc( 1 , fltarr ( 11 ) ) 
si it  =  sl ( 1 ) 

vk  =  assoc( 3,fltarr( 512)  ) 
wl=indgen( 512 ) *3 . 133  +1798 
i  =  0 

for  t=100, 2000, 100  do  begin 
i  =  i  +  l 


int=wl*0 
intl=wl *0 
int2=wl *0 
int3=wl*0 
in 1 4=wl *  0 
i nt 5=wl *  0 
i nt6=wl *  0 
int7=wl*0 
int8=wl*0 
int9=wl*0 
intl0=wl*0 
intll=wl *0 
intl2=wl*0 
inti 3=wl *  0 
intl4=wl *0 
intl5=wl*0 

i nt 1 6=wl *  0 
i nt 1 7  =wl *  0 
intl8=wl *0 
intl9=wl*0 
int20=wl *0 
int21=wl *0 

int22=wl*0 

int23=wl*0 
int24=wl*0 
int25=wl *0 
i nt26=wl *  0 
int27=wl*0 
i n 1 2  8=wl *  0 


Calculates  synthetic  Vegard-Kaplan 
spectra  at  temps  from  100  to  2000 
degrees  for  wavelengths  from  1800  to 
3400  A.  Results  convolved  with  MUSTANG 
instrument  slit  function 


;  these  variables  provide  a  means  of  tweaking 
;  individual  transitions 


t-  calculates  the  transitions'  band  systems 

into  the  total  spectrum. 


TRANSITION 


WL 


summing  each  one 


vksyn , t , 6 , 1 , wl , int  ;  1799 
vksyn , t , 4 , 0 , wl , int  ;  1808 
vksyn , t , 1 3 , 5 , wl , int  ;  1823 
vksyn , t , 1 1 , 4 , wl , i n t  ;  1824 
vksyn, t, 9, 3, wl, int  ;  1827 
vksyn , t , 7 , 2 , wl , int  ;  1833 
vksyn,t,5,l,wl,int  •  1841 
vksyn , t , 3 , 0 , wl , int  ;  1853 
vksyn, t, 12, 5, wl, int  ;  i860 
vksyn, t, 10, 4 ,wl , int  ;  18^3 
vksyn, t, 8, 3, wl, int  ;  1868 
vksyn, t, 6, 2, wl, int  ;  1876 
vksyn , t , 4 , 1 , wl , int  •  1887 
vksyn , t , 1 3 , 6 , wl , int  •  1899 
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vksyn, t , 11 , 5 , wl , int  ;  1900 
vksyn,t,2,0,wl,int  ;  I901 
vksyn, t, 9 , 4 ,wl , int  ;  1905 
vksyn, t , 7 , 3 , wl , int  ;  1912 
vksyn, t, 5,2, wl, int  ;  1923 
vksyn, t, 3 , 1 ,wl , int  ;  1935 
vksyn, t , 12 , 6 ,wl , int  ;  1939 
vksyn, t , 10 , 5 , wl , int  ;  1943 
vksyn ,t,8,4,wl,int  ;1950 
vksyn , t , 1 , 0 ,wl , inti  ;  1954 
int=int+intl*l .  5 

vksyn, t , 6 , 3 ,wl , int  ;  i960 
vksyn , t , 4 , 2 , wl , int  ;  1973 
vksyn, t , 13 , 7 ,wl , int  ;  1980 
vksyn , t , 11 , 6 , wl , int  ;  1983 
vksyn,t,9,5,wl,int  ;  1989 
vksyn, t , 2 , 1 ,wl , int  ;  1990 
vksyn  ..  t ,  7 , 4  ,  wl ,  int  ;  1998 
vksyn, t,0,0,wl,int2  ;  2OIO 
int=int+int2*2.6 

vksyn, t, 5, 3, wl, int  ;  2OII 
vksyn, t, 12, 7, wl, int  ;  2024 
vksyn, t , 3 , 2 , wl , int  ;  2027 
vksyn, t, 10, 6, wl, int  ;  2029 
vksyn , t , 8 , 5 , wl , int  ;  2038 
vksyn , t , 1 , 1 ,wl , int3  ;  2047 
int=int+int3*l . 5 

vksyn, t , 6 , 4 , wl , int  •  2050 
vksyn, t, 4 , 3 ,wl , int  ;  2065 
vksyn , t , 13 , 8 , wl , int  ■  2067 
vksyn, t , 11 , 7 , wl , int  ;  2071 
vksyn, t , 9 , 6 , wl , int  ;  2079 
vksyn , t , 2 , 2 , wl , int  ;  2085 
vksyn, t , 7 , 5 , wl , int  ;  2090 
vksyn, t, 5, 4, wl, int  ;  2106 
vksyn, t,0,l,wl,int4  ;  2109 
int=int+int4*2.6 

vksyn,t,12,8,wl,int  ;  2115 
vksyn, t, 10, 7, wl, int  ;  2122 
vksyn , t , 3 , 3 , wl , int  ;  2125 
vksyn, t, 8 , 6 ,wl , int  ;  2133 
vksyn , t , 6 , 5 , wl , int  ;  2147 
vksyn, t, 1 , 2 ,wl , int5  ;  2148 
int  =  int+int5*l .  5 

vksyn , t , 13 , 9 , wl , int  ;  2161 
vksyn , t , 4 , 4 , wl , int  ;  2166 
vksyn, t , 11 , 8 , wl , int  ;  2167 
vksyn, t , 9 , 7 ,wl , int  ;  2177 
vksyn , t , 2 , 3 , wl , int  ;  2189 
vksyn , t , 7 , 6 , wl , int  ;  2191 
vksyn , t , 5 , 5 , wl , int  ;  2209 
vksyn, t , 12 , 9 ,wl , int  ;  2213 
vksyn , t , 0 , 2 ,wl , int6  ;  2216 
int=int+int6*2.6 

vksyn, t, 10, 8, wl, int  ;  2222 
vksyn , t , 3 , 4 , wl , int  ;  2231 
vksyn , t , 8 , 7 , wl , int  ;  2236 
vksyn, t, 6, 6, wl, int  ;  2253 
vksyn , t , 1 , 3 , wl , int7  ;  2258 
int  =  int  +  int7  *1.5 
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vksyn,t,ir,10,wl,int 

;  2262 

vksyn,t,ll,9,wl,int 

;  2270 

vksyn , t , 4 , 5 , wl , int 

;  2275 

vksyn , t , 9 , 8 , wl , int 

;  2282 

vksyn, t , 7 , 7 ,wl , int 

;  2299 

vksyn, t , 2 , 4 ,wl , int 

;  2302 

vksyn, t , 1 2 , 1 0 , wl , i nt 

;  2320 

vksyn, t , 5 , 6 ,wl , int 

;  2321 

vksyn, t, 10,9, wl , int 

;  2331 

vksyn, t , 0 , 3 , wl , int8 

;  2334 

int=int+int8*2.6 

vksyn, t , 3 , 5, wl , int 

;  2347 

vksyn , t , 8 , 8 , wl , int 

;  2347 

vksyn, t, 6 , 7 ,wl , int 

;  2368 

vksyn, t , 13 , 11 , wl , int 

;  2372 

vksyn, t , 1 , 4 , wl , int9 

;  2379 

int=int+int9*1.5 

vksyn, t, 11, 10, wl, int 

2382 

vksyn, t , 4 , 6 , wl , int 

2394 

vksyn, t, 9, 9, wl, int 

2397 

vksyn, t , 7 , 8 , wl , int 

2417 

vksyn, t, 2 , 5 ,wl , int 

2425 

vksyn, t , 12 , 11 ,wl , int 

2435 

vksyn, t, 5, 7, wl, int 

2443 

vksyn, t , 10 , 10 , wl , int 

2449 

vksyn,t,0,4,wl,intl0 

2463 

int  =  int  +  intl0*2 . 6 

vksyn, t , 8 , 9 , wl , int 

2469 

vksyn,  t ,  'i ,  6  ,  wl ,  int 

2474 

vksyn, t, 13 , 12 ,wl , int 

2490 

vksyn, t, 6, 8, wl, int 

2494 

vksyn, t, 11, 11, wl, int 

2504 

vksyn, t , 1 , 5 ,wl , intll 

2511 

'  int  =  int  +  intll*l .  5 

vksyn, t, 9, 10, wl, int 

2522 

vksyn , t , 4 , 7 , wl , i nt 

2524 

vksyn, t , 7 , 9 , wl , int 

2547 

vksyn , t , 2 , 6 , wl , int 

2561 

vksyn, t , 12 , 12 , wl , int 

2561 

vksyn , t , 5 , 8 , wl , int 

2577 

vksyn, t, 10 , 11 ,wl , int 

2578 

vksyn, t, 8 , 10 ,wl , int 

2602 

vksyn , t , 0 , 5 , wl , intl2 

2605 

int=int+inti2*2.6 

vksyn , t , 3 , 7 , wl , int 

2613 

vksyn, t , 13 , 13 ,wl , int 

2620 

vksyn , t , 6 , 9 , wl , intl3 

2632 

int=int+intl3* . 83 

vksyn, t, 11 , 12 ,wl ,  int 

2637 

vksyn , t , 1 , 6 , wl , intl4 

2657 

int=int+intl4*i,5 

vksyn, t , 9 , 11 , wl , int 

2659 

vksyn , t , 4 , 8 , wl , intis 

2668 

int=int+intl5*.67 

vksyn, t , 7 , 10 ,wl , intl6 

2689 

int=int+intl6* . 6 

vksyn, t , 12 , 13 ,wl , int 

2698 

vksyn , t , 2 , 7 , wl , int  ; 

2711 

vksyn, t , 10 , 12 ,wl , int 

2719 

vksyn, t, 5, 9, wl, int 

2724 
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vksyn, t,8,ll,wl,int 

;  2748 

vksyn,t,13,14,wl,int 

;  2761 

vksyn ,t,0,6,wl,intl7 

;  2762 

int=int+intl7*2.6 

vksyn , t , 3 , 8 , wl , int 

;  2767 

vksyn,t,ll,13,wl,intl8 

;  2782 

int=int+intl8* . 7  1 

vksyn, t , 6 , 10 , wl , int 

;  278  3 

vksyn,t,9,12,wl,int 

;  2810 

• 

vksyn,t,l,7,wl,intl9 

;  2818 

int  =  int  +  intl9*l .  5 

vksyn , t , 4 , 9 , wl , int 

;  2826 

vksyn, t, 7 , 11 ,wl , int 

;  2845 

- 

vksyn ,t,12,14,wl,int20 

;  28  4  8 

int=int+int20* ,67 

vksyn,t,10,13,wl,int 

;  2875 

vksyn,t,2,8,wl,int21 

;  2877 

int=int+int21* . 67 

vksyn,t,5,10,wl,int22 

:  2887 

int=int+int22* . 5  1 

vksyn, t,8,12,wl,int23 

;  2909 

( 

1 

int=int+int23* . 5 

vksyn ,t,13,15,wl,int 

;  2916 

vksyn ,t,0,7,wl,int24 

:  2937 

int  =  int  +  int24*2 . 6 

vksyn , t , 3 , 9 , wl , int 

2938 

vksyn, t , 11 , 14 , wl , int 

2942 

vksyn, t, 6 , 11 ,wl , int 

2951 

vksyn , t , 9 , 1 3 , wl , int 

2976 

vksyn, t, 1 , 8 ,wl , int25 

2998 

int  =  int  +  int25*l .  5 

vksyn, t , 4 , 1 0 , wl , i nt 

3002 

vksyn,t,12,15,wi,int 

3013 

vksyn ,t,7,12,wl,int 

3017 

% 

vksyn , t , 10 , 14 , wl , int 

3046 

vksyn , t , 2 , 9 , wl , int 

3062 

I 

vksyn , t , 5 , 1 1 , wl , i nt 

3068 

vksyn, t, 13, 16, wl, int 

3087 

vksyn , t , 8 , 1 3 , wl , int 

3087 

vksyn, t, 11 , 15 ,wl , int 

3119 

vksyn , t , 3 , 1 0 , wl , i nt 

3128 

vksyn,t,0,8,wl,int26 

3133 

int=int+int26*2.6 

vksyn , t , 6 , 1 2 , wl ,  i nt 

3137 

vksyn , t , 9 , 14 , wl , int 

3159 

vksyn, t , 12 , 16 ,wl , int 

3196 

vksyn , t , 4 , 11 , wl , int 

3198 

vksyn ,t,l,9,wl,int27 

3200 

int=int+int27*1.5 

vksyn , t , 7 , 1 3 , wl , int 

3210 

vksyn , t , 1 0 , 1 5 , wl , i nt 

3235 

vksyn , t , 2 , 10 , wl , int 

3269 

vksyn , t , 5 , 12 , wl , int 
vksyn,t,13,17,wl,int 

3270 

3276 

vksyn , t , 8 , 1 4 , wl , int 

3285 

vksyn , t , 1 1 , 1 6 , wl , int 

3315 

vksyn , t , 3 , 11 ,wl , int 

3342 

vksyn , t , 6 , 13 , wl , int 

3345 

vksyn,t,0,9,wl,int28 

3353 

- 

int=int+int28*2 . 6 
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vksyn , t , 9 , 1 = . wl , ' nt 
vksyn , t , 12 , 17 , wl , int 

vk( i )=convol( int,slit) 


endf o  r 
do?  ? ,  1 
closed 
return 
end 

@ [ synspec ] vksyn 


;  3  36  4 
;  3399 
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;Walden, 1991 

PRO  NODLTCALC,WL,INT,SLIT 


OPENU ,10,' NODLT . DAT ' 
DLT=ASSOC( 10 , FLTARR( 512 ) ) 

WL=INDGEN(  512  ) ’"3 . 133  +  1798 


;This  procedure  calculates  synthetic 
;  nitric  oxide  delta  band  spectra 
;  at  temps  from  100  to  2000  deg  and 
;  convolves  results  with  MUSTANG  instru- 
;  ment  slit  function. 


OPENR, 11 , ' [WALDENJSLIT.DAT' 

SL=ASSOC( 11 , FLTARR( 11 ) ) 

SLIT=SL( 1 ) 

1  =  0 

FOR  T=100 , 2000 , 100  DO  BEGIN 

1=1+1  ;  Wavelength 


INT=WL*0 

DELSyN,T,0,0,WL,INT  ;  1909.8 

DELSyN,T,0,l,WL,INT  ;  1980.5 

DELSyN,T, 0 , 2 ,WL, INT  ;  2055.9 

DELSyN,T,0,3,WL,INT  ;  2134.8 

DELSyN.T, 0 , 4 ,WL, INT  ;  2219.6 

DELSyN,T, 0 , 5 ,WL, INT  ;  2310.1 

DELSyN,T,0,6,WL,INT  ;  2407.6 

DELSyN,T,0,7,WL,INT  ;  2509 

DELSyN,T,0,8,WL,INT  ;  2620 

DELSyN,T,0,9,WL,INT  ;  2739 

DELSyN,T,0,10,WL,INT  ;  2866 

DELSyN,T, 0 , 11 ,WL, INT  ;  3003 

DELSyN,T,0,12,WL,INT  ;  3152 

DELSyN, T, 0 , 13 ,WL, INT  ;  3313 

DELSyN,T, 0 , 14 ,WL, INT  ;  3488 


DLT( I )=CONVOL( INT, SLIT) 

ENDFOR 

CLOSE, 10 
CLOSE, 11 

'.ETURN 

;nd 

[ SyNSPEC]DELSyN 
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;Walden ,1991 

;This  program  derives  the  synthetic  spectrum  using  tweaking  factors 
, -first  determined  by  Mack(1991)  for  1800-2100A,  and  extended  the  full  range 
;of  the  MUSTANG  data. 

PRO  NOGAMCALC,WL, INT, SLIT 


OPENU ,10,' NOGAM . DAT ' 

GAM=AS  SOC (10, FLTARR (512)) 
WL=INDGEN( 512  )  "3 . 133  +  1798 


;This  procedure  calculates  synthetic 
;  nitric  oxide  gamma  band  spectra 
;  at  temps  from  100  to  2000  deg  and 
;  convolves  results  with  MUSTANG  instru- 
;  ment  slit  function. 


OPENR, 11 , ' [WALDEN]SLIT.DAT' 

SL=ASSOC( 11 ,FLTARR( 11 ) ) 

SLIT=SL(  1 ) 

1  =  0 

FOR  T=100 , 2000 , 100  DO  BEGIN 
1  =  1  +  1 

;These  variables  allow  for  the  adjusting  of 

;individual  transition's  probability 

INT=WL*0 

INT1=WL*0 

INT2=WL*0 

INT3=WL*0 

INT4=WL*0 

INT5=WL*0 

INT6=WL*0 

INT7=WL*0 

INT8=WL*0 

INT9=WL*0 

INT10=WL*0 . 

INT11=WL*0  . 

INT12=WL*0  . 

INT13=WL*0 . 

INT14=WL*0 . 

INT15=WL*0 . 

INT16=WL*0 . 

INT17=WL*0 . 

INT18=WL*0  . 

INT19=WL*0 . 

INT20=WL*0 . 

INT21=WL*0 . 

INT22=WL*0 . 

INT23=WL*0 . 

This  section  calculates  the  spectrum  from  each  transition's  band 

system . 

Transition  WL 


GAMSYN,T,3,0,WL,INT 
GAMSYN, T, 3 , 1 ,WL, INT 
GAMSYN,T, 2 , 0 ,WL, INT9 
INT=INT+INT9* .75 
GAMSYN, T, 3 , 2 ,WL, INT 
GAMSYN,T,2,1,WL,INT 
GAMSYN, T,1 ,0,WL, INTI 
INT=INT+INT1+1 .125 
GAMSYN, T , 3 , 3 , WL , INT 
GAHSYN,T,2,2,WL,INT2 


1956 . 1 

2029.8 

2046.5 

2108 . 5 

2128.4 

2148.1 

2192.8 

2215.4 
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INT=INT+INT2*1 . 26 
GAMSYN,T, 1 , 1 ,WL, INT3 
INT=INT+INT3*1 . 3 
GAMSYN, T, 0 , 0 ,WL, INT4 
INT=INT  +  INT4*1 . 2 
GAMSYN, T, 3 , 4 ,WL, INT 
GAMSYN, T, 2 , 3 ,WL, INT5 
INT-INT+INT5*! . 6 
GAMSYN, T, 1 , 2 ,WL, INTIO 
INT=INT+INT10*1 . 8 
GAMSYN, T,0,1,WL,INT6 
INT=INT+INT6*1 . 2 
GAMSYN, T, 3, 5, WL, INT 
GAMSYN, T, 2, 4, WL, INT 
GAMSYN, T,l, 3, WL,INT7 
INT=INT+INT7*1 . 6 
GAMSYN , T , 0 , 2 , WL , INT8 
INT=INT+INT8*1 . 2 
GAMSYN,T,3,6,WL,INT 
GAMSYN, T, 2, 5, WL, INT 
GAMSYN, T, 1 , 4 ,WL , INTI  7 
INT=INT+INT17*1 . 3 
GAMSYN, T,0, 3, WL, INTll 
INT=INT+INT11*1 . 2 
GAMSYN, T, 3 , 7 ,WL, INT 
GAMSYN, T, 2, 6, WL, INT 
GAMSYN, T, 1 , 5 ,WL , INTI 2 
INT=INT+INT12*1 . 3 
GAMSYN, T, 3, 8, WL, INT 
GAMSYN , T , 0 , 4 , WL , INTI  3 
INT=INT+INT13*1 . 2 
GAMSYN, T, 2, 7, WL, INT 
GAMSYN, T, 1 , 6 ,WL , INTI 4 
INT=INT+INT14*1 . 3 
GAMSYN, T, 3, 9, WL, INT 
GAMSYN, T,0, 5, WL,INT15 
INT=INT+INT15*1 . 2 
GAMSYN, T, 2, 8, WL, INT 
GAMSYN, T, 1 , 7 ,WL, INT16 
INT=INT+INT16*1 . 3 
GAMSYN, T, 3, 10, WL, INT 
GAMSYN, T,0, 6, WL,INT18 
INT=INT+INT18*1 . 2 
GAMSYN, T, 2, 9, WL, INT 
GAMSYN,T, 1 ,8 ,WL, INT19 
INT=INT+INT19*1 . 3 
GAMSYN, T, 3, 11, WL, INT 
GAMSYN,T, 0 ,7 ,WL, INT20 
INT=INT+INT20*1 . 2 
GAMSYN , T , 2 , 1 0 , WL , I NT 
GAMSYN,T, 1 ,9 ,WL, INT21 
INT=INT+INT21*1 , 3 
GAMSYN, T, 3, 12, WL, INT 
GAMSYN, T, 2, 11, WL, INT 
GAMSYN, T, 0 , 8 ,WL, INT22 
INT=INT+INT22*1 . 2 
GANSYN,T,3,13,WL,INT 
GAMSYN, T,l, 10, WL,INT23 
INT=INT+INT23*1 , 3 


2238 . 3 
2261.8 

2283.0 

2308 . 5 

2335 

2362.2 

2378.8 
2409 

2438.6 

2470.1 

2481.2 

2516.4 

2550.3 

2586.3 

2591 

2629.9 

2670.1 

2709 

2712.0 

2753.6 
2800.0 

2836 

2848.2 

2888.2 

2940.5 

2973 

2997.6 

3044.3 

3112.4 

3120.6 

3170.7 

3201.1 

3278.5 

3303.0 

3361 

3375.5 

3456 

3458 . 5 
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GAM( I )=CONVOL( INT, SLIT) 


ENDFOR 

CLOSE, 10 
CLOSE, 11 

RETURN 

END 

@[SYNSPEC]GAMSYN 


;Walden , 1991 
PRO  NOEPSCALC 


OPENU , 1 , ' NOEPS . DAT ' 
EPS=ASSOC( 1 , FLTARR( 512 ) ) 


WL=INDGEN( 512 ) *3 . 133  +  1798 


This  procedure  calculates  synthetic 
nitric  oxide  epsilon  band  spectra 
at  temps  from  100  to  2000  deg  and 
convolves  results  with  MUSTANG  instr- 
ment  slit  function,  where  EPS{1) 
corresponds  to  T=100,etc.  It  uses 
scale  factors  determined  by  Bosserman, 
except  for  (0,1),  which  Bosserman 
used  .84  to  tweak 


OPENR,2, ' [WALDEN1SLIT.DAT' 
SL=ASSOC( 2 , FLTARR( 11 ) ) 
SLIT=SL( 1 ) 


1  =  0 

FOR  T=100 , 2000 , 100  DO  BEGIN 
1  =  1  +  1 
INT=WL*0 
INT1=WL*0 
INT2=WL*0 
INT3=WL*0 
INT4=WL*0 

EPSSYN,T,1,0,WL,INT 

EPSSYN,T,2,2,WL,INT 

EPSSYN,T,1,1,WL,INT 

EPSSYN,T,0,0,WL,INT1 

INT=INT+INT1* . 63 

EPSSYN,T,2,3,WL,INT 

EPSSYN,T,1,2,WL,INT 

EPSSYN,T,0,1,WL,INT2 

INT=INT+INT2* .90 

EPSSYN,T,1,3,WL,INT 

EPSSYN,T,0,2,WL,INT3 

INT=INT+INT3*1 .03 

EPSSYN,T,1,4,WL,INT 

EPSSYN,T,0,3,WL,INT4 

INT=inT+INT4*l .11 

EPSSYN,T,1,5,WL,INT 

EPSSYN,T, 0 , 4 ,WL, INT 

EPSSYN,T,2,8,WL,INT 

EPSSYN,T,2,9,WL,INT 


;  Wavelength 


;  1799.6 
;  1849.5 
;  1863.2 
;  1876.6 

;  1918 
;  1932.3 
;  1945.0 

;  1998.6 
;  2017.5 

;  2073.0 
;  2  09  4.5 

;  2157.5 
;  2176.1 
;  2294.0 
;  2385.0 


EPS( I )=CONVOL( INT, SLIT) 

ENDFOR 

CLOSE, 2 
CLOSE, 1 

ETURN 

ND 

[ synspec ] EPSSYN 
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;Walden ,1991 
pro  01  2972calc 


;This  program  creates  a  synthetic  model  for  a 
;  neutral  atomic  oxygen  emission 


openu,l, 'OI_2972.dat' 
dat=assoc( l,fltart{512)) 


openr , 2 , ' (waldenjslit.dat' 
sl  =  assoc(2,fltarr(ll)  ) 

slit-sl(4)  ;  adjusts  wl  to  2972  vice 


2972.875 


hold=f Itar r ( 512 ) 


hold=indgen(512)*0.0 

hold(375)=1.0 


dat(0).convol(hold,slit)/3.133  ; spec t ra-l/wl  bin  so  that  scale  factor  is 

/•absolute  intensity 

close  ,  2 
close ,  1 


return 

end 
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;program  calculates  synthetic  spectrum 
;for  0+  atomic  emission  @  2470A 


;Walden ,1991 
pro  OII_2470calc 

openu,l,'0II  2470.dat' 
dat=assoc(l,Tltarr(512)  ) 

openr,2, ' (waldenlslit.dat' 
sl=assoc(2,fltarr(ll)  ) 
slit=sl(l)  ;this  makes  wl  centered  @  2968. 5A 

hold=f Itar r ( 512  ) 

hold=indgen(512)*0.0 
hold( 214 )=1 .0 

da t ( 0 ) =convol ( hold , si i t )/3 , 1 33  ;spectra  is  1/wl  bin  so  scale  factor 

;is  absolute  intensity 


close ,  2 
close ,  1 

return 

end 
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;Walden ,1991 

;This  program  creates  a  syntheic  model  for  the  0++  atomic  emission  @  2853A. 
pro  OIII_2853calc 

openu , 1 , ' 01 1 1  2853  .dat' 
dat=assoc(l,fTtarr(512) ) 

openr , 2 , ' [ walden ] si i t . dat ' 

.  sl  =  assoc ( 2 , f Itar r  ( 11 ) ) 
si i t  =  sl ( 1 ) 

hold=f Itar r ( 512  ) 

hold=indgen(512)*0.0 
hold( 337 )=1 .0 

spike=convol(hold,slit) 

dat ( 0 )=spike/total ( spike )/3 . 133  ;dat(0)  is  regular  convolved  spike 

spike ( 338 : 510 )=spike ( 339 ; 511 )  ;dat(l)  is  spike  with  adjustment  for  apparent 

;  dropped  point  in  data. 
dat( 1 )=(spike/total(spike ) )/3.133 

close ,  2 
close  ,  1 

return 

3nd 
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;Walden ,1991 

;This  program  calculates  a  synthetic  spectrum  for  the  Nil  doublet 
;with  the  ratio  based  on  the  findings  of  Bucsela  and  Sharp  (1989) 
pro  NlI_doubcalc 

openu,l,'NII  2143doub.dat' 
nii  =  assoc(l,lltarr(512)  ) 

ope nr, 2, ' [waldenlslit.dat' 
sl  =  assoc(2,fltarr(ll)  ) 

hold=fltarr( 512) 

hold( 109)=. 58  ;2139,68A  line 

intl  =  convol ( hold , si ( 1 )  )  ;wl  is  in  center  of  pixel 

hold2=fltarr( 512) 

hold2(110)=1.0  ;2143.55  line 

int2=convol ( hold2 , si ( 3 ) ) ; wl  is  1/4  to  the  right  of  center  of  pixel 
int=intl+int2 

nii(l)=(int/total(int) )/3.133 

;total  of  spectral  area  is  1/wl  bin,  so 
;scale  factor  will  be  the  absolute  intensity 


close ,  2 
close ,  1 
return 
end 
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;V3alden,1991 

;This  program  creates  a  synthetic  spectrum  of  atomic  emissions  observed  in 
;the  March  1990  MUSTANG  data, 
pro  otherspec 

ope nr, 1, ' [walden]slit.dat' 
slit=assoc(l,fltarr(ll) ) 

•  openu , 2 , ' othe r spec . da t ' 
spect=assoc(2,fltarr(512) ) 

spec  =  f 1  tar  r ( 51 2  ) 

■  numlines=20  jnumber  of  contributing  atomics 

;wavelengths  determined  from  Bureau  of  Standards 
, -listing  of  atomic  emissions  and  NASA  listing  of  atomic 
; emi s  s i ons 


Wavelength 

Species 

Pixel 

1846 . 0 

NIII 

15 

2182 . 6 

on 

123 

2190.6 

on 

125 

2337.8 

ARI I 

172 

2384.6 

OIV 

187 

2445 . 6 

on 

207 

2490 . 3 

NI I 

221 

2509 . 9 

Nil 

227 

2665.8 

OIII 

277 

2709.8 

Nil 

291 

2741  .  3 

on 

301 

2885 . 3 

NI  I 

347 

2973.6,2974.6, 

Nil  multiplet 

375-377 

2975.7,2975.9, 

2976.9,2977.3, 

2978.4,2978.6 


nx= [15, 123, 125, 172, 187, 207, 221, 227, 277, 291, 301, 347,$ 
375,375,376,376,376,376,377,377] 

nt=(10, 2, 1,10, 9, 3, 3, 4, 5, 6, 4, 4, 2, 3, 3, 4, 4, 4, 5, 5] 

>lno=[ 3, 4, 3, 3, 4, 4, 1,3, 1,1, 1,1, 3, 2, 4, 1,3, 2, 4, 4] 

or  i=0,numlines-l  do  begin 
blank=fltarr(512) 
blank(pix(i) )=.l*int(i) 
impulse=convol(blank,slit(slno( i ) ) ) 
spec=spec+ impul se 

ndf  or 

pect( 0)=indgen( 512 1*3.133+1798 

pec t ( 1 ) = ( spec/tota 1 { spec ) ) /3 . 1 3 3  ;so  that  scale  factor  is  absolute  int. 


lose  ,  2 
lose  ,  1 
e  tu  r  n 
nd 


9.5 


APPENDIX  C 


Appendix  C  is  the  programs  used  to  determine  the  standard 
deviations  of  the  means (sdoms)  for  the  data.  The  sdoms  are 
used  to  determine  the  weighting  factors  for  the  fitting 
routine  by  the  program  titled  "weights2,"  which  is  also 
included . 
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; Walden ,1991 

pro  m90sdev , al t , indexl , index2 

;This  procedure  calculates  the  std.dev  and  std  dev  of  the  mean 
;for  the  averaged  spectra. 

;Inputs  consist  of  the  requested  altitude  (1-25)  and  its  corresponding 
, -beginning  and  ending  index  numbers. 

/Outputs  are  the  SDOM  (unshifted  and  shifted)  and  are  written  to 
; um9  0_sdom . dat  and  m90_sdom.dat,  respectively 


navg=0 
st= '  ' 

sigma=dblarr(512) 
sigmam=dblarr( 512  ) 

openr,l, 'dubO: [claytonlmustang_raw.dat' 
dat  =  assoc ( 1 , intar r ( 512  )  )  ~ 

openr,2, ' [mustang] um9  0_avg . da t ' 
avg=assoc(2,fitarr(512)  ) 

openu,3,  '  [mustangjum9  0_sdrri.dat' 
sdom=assoc(3,fltarr(512)  ) 


/opens  the  raw  data  file 

/ opens  the  averaged  spectrum  file 
/  (unshifted  version) 

/opens  the  file  to  write  SDOM  into 


T=index2-indexl 
for  j=0,n  do  begin 

print, indexl+ j 

print, 'Do  you  want  to  include  this  record' 
print, 'Type  <CR>  for  yes,  N  for  no.' 
read , st 

if  St  eq  'n'  then  goto,labell 
'  sigma=sigma  +  ( da t ( i ndexl+ j ) -avg ( al t ) ) ' 2 

navg=navg+l 
labell : dummy=0 

?ndf or 
)ree=navg-l 

iigma=sqrt( sigma/ free) 

,igmam=sigma/sqrt(navg) 

.dom( alt ) =sigmam 

;  1  o  s  e  ,  1 
: lose  ,  2 
lose , 3 

eturn 

nd 
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; Walden ,1991 
pro  weights2 

;This  program  checks  for  values  where  the  detector  saturated,  and  gives 
;weights  based  on  the  deviations  of  the  mean  being  equal  to  the  data  at 
, 'these  points.  It  also  gives  a  small  weight  factor  to  those  elements 
/corresponding  to  nitrogen  second  positive  wavelengths  (sets  sdom  equal  to 
/largest  data  value  in  N2  2nd  pos  bands).  Since  sdom  ^as  given  value  of  data  at 
/dropped  points,  it  gives  small  weights  to  dropped  values  (with  the  exception 
/of  2972A) 

openr,l, ' [mustanglm90  intsdm.dat' 
dom=assoc (l,fltarr(51^) ) 

ope nr , 2 , ' [ mustang ] m90_int . da t ' 
dat  =  assoc(2,fltarr(512)  ) 

openu , 3 , ' m90_domwt . dat ' 
domwt=assoc ( 3 , f Itar  (512)) 

openr,4,'[mack] mack_cal . da t ' 
ss  =  assoc(4,fltarr(512)  ) 
sens  =  ss ( 22  ) 

ope nr, 5, '  [mustang] m90_avg . dat ' 
dd  =  assoc(5,fltarr(512)  ) 

big=1800  /approx  N2  2nd  pos  largest  datapoint  value 

holddom=f Itar r ( 512  ) 

iolddomwt=fltarr(5i2) 

•iolddat  =  fltarr(  512) 

3at=900 . 00/sens 

for  i=l,25  do  begin 

holddom=dom ( i  ) 
holddat=dat ( i  ) 
for  j=0,511  do  begin 

; i f  detector  was  saturated,  make  sdom  equal  to  data  value 

if  (holddat(j)  gt  sat(j))  then  holddcm ( j ) =holdda t ( j ) 


if 

wavelength 

i  s 

in  N2 

2nd 

pos 

t 

make 

sdom  equal  to  big 

wavelength 

f 

(  j 

eq 

277  ) 

or 

(  j 

eq 

278  ) 

or 

(  j 

eq 

279) 

then 

holddem( j ) =big 

2669 

f 

(  j 

eq 

318  ) 

or 

(  j 

eq 

319  ) 

or 

(  j 

eq 

320  ) 

then 

holddom( j )=big 

2796 

f 

(  j 

eq 

322  ) 

or 

(  j 

eq 

323  ) 

or 

(  j 

eq 

324  ) 

then 

holddom( j )=big 

2811 

f 

(  j 

eq 

362) 

or 

(  j 

eq 

363) 

or 

(  j 

eq 

364) 

then 

holddom( j ) =bi g 

2934 

f 

(  j 

eq 

368  ) 

or 

(  j 

eq 

369  ) 

or 

(  j 

eq 

370) 

then 

holddom( j ) =big 

2953 

f 

(  j 

eq 

374  ) 

or 

(  j 

eq 

375) 

or 

(: 

eq 

376) 

then 

holddom( j )=big 

2973 

f 

(  j 

eq 

409  ) 

or 

(  j 

eq 

410) 

or 

( j 

eq 

411 ) 

then 

holddom( j ) =big 

3083 

f 

(  j 

eq 

417  ) 

0  r 

(j 

eq 

418  ) 

or 

( j 

eq 

419  ) 

then 

holddom( j ) =bi g 

3107 

f 

(  j 

eq 

425) 

or 

(  j 

eq 

426  ) 

or 

( j 

eq 

427  ) 

then 

holddom( j ) =big 

3132 

f 

(  j 

eq 

433  ) 

or 

(  j 

eq 

434  ) 

or 

( j 

eq 

435) 

then 

holddom( j ) =big 

3157 

f 

(  j 

eq 

461  ) 

or 

(  j 

eq 

^62) 

or 

( j 

eq 

463) 

then 

holddom( j )=big 

3244 

f 

(  j 

eq 

470  ) 

o  r 

(  j 

eq 

471  ) 

or 

( j 

eq 

472  ) 

then 

holddom( j )=big 

3274 

f 

(  j 

eq 

480  ) 

or 

(  j 

eq 

481  ) 

or 

( j 

eq 

482  ) 

then 

holddom( j ) =bi g 

3305 

f 

(  j 

eq 

490  ) 

o  r 

(  j 

eq 

491  ) 

0  r 

( j 

eq 

492) 

then 

holddom( j ) =big 

3337 

f 

(  j 

eq 

501  ) 

or 

(  j 

eq 

502  ) 

or 

( j 

eq 

503) 

then 

holddom( j ) =big 

3370 

f 

(  j 

gh 

493  ) 

then 

holddom( j )=big 

/last  spike 

now  give  weights  their  values 

if  holddom(j)  eq  0.  then  holddomwt ( j ) =1 . 
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if  holddom(j)  ne  0.  then  holddomwL ( j ) =1 . /holddom ( j ) " 2 

endf or 

/Correct  for  where  low  value  of  data  point  at  dropped  word  gives  false 
/high  weight  to  point  by  giving  dropped  points  1/10  the  maximum  weight 
/of  any  other  point 

/(except  for  2972A  -pixels  371  to  376  and  2850A  -333  to  336). 
dropped=whe re ( dd ( i )  eq  -1.) 
hold=holddomwt 
hold ( dropped ) =0 . 0 
holddomwt (dropped)=.l*max(hold) 

/  holddomwt ( 37 1 : 376  ) =max ( hold  ) 

/  holddomwt ( 33 3 : 3 36 ) =max ( hold ) 

/assign  weight  to  file  m90_domwt.dat 
domwt ( i ) =holddomwt 

endf or 
close ,  1 
close  ,  2 
close ,  3 
close ,  4 
close ,  5 

return 

end 
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a  ij 


;Walden , 1991 

;This  program  shifts  the  MUSTANG  data  array  according  to  dropped  pixels  as 
;input  in  the  procedure  shiftarray.  ^  ' 

pro  shii.i.data 

openr , 1 , ’ [ mustang ] um9 0_avg . dat ' 
mm=assoc (l,fltarr(512)) 

openu , 2 , ' [ mustang ]m90_avg .dat ' 
dd=assoc (2,fltarr(512)  ) 

openr , 3 , ' [mustanglum90_sdm.dat' 
uu=assoc (3,fltarr(512)) 

openu , 4 , ' [ mustanglm90_sdm.dat' 
ss=assoc (4,fltarr(512) ) 

wl=indgen( 512) *3. 133+1798 


dd( 0 )=wl 
ss ( 0 ) =wl 

for  i=l,25  do  begin 

shiftarray, mm (i),holdl 
dd( i )=holdl 

shiftarray, uu{i) ,hold2 
ss ( i ) =hold2 

endf or 

lose ,  4 
lose ,  3 
:lose ,  2 
:lose ,  1 

■return 

?nd 
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; Walden ,1991 

;This  program  shifts  an  array  of  data  according  to  the  "drop"  array  of 
;  p  i  X  e  1  s  . 

pro  shi f tar  ray , ina r r , outar r 

dr op=[ 46 2, 4 4 6, 4 3 0,4 14, 3 98, 382, 3 66, 3 5 0,3 50, 3 34, 3 18, 3 02, 2 8 6, 27 0,2 5 4, 2 3 9, 2 2 3, 22  3  ,$ 
207,191,175,159,143,127,111,95,95,80,80,64,48,32,16] 

ndrops=n_e laments ( drop ) 
outar r=inarr 

length=n_elements ( inar  r ) 
for  i=0,ndrops-l  do  begin 

;  This  section  is  to  put  a  -1  in  dropped  values  and  shift  as  required. 

outarr{drop(i)+l)=outarr(drop(i):length-2) 
outarr(drop( i )  )=-l  .0 

endf or 

return 

end 
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; Walden , 1991 

;This  program  takes  the  MUSTANG  data  from  the  averaged,  shifted 

;and  puts  in  the  appropriate  average,  and  adjusts  for  instrument 
pro  average 


records 
response . 


openr , 1 , ' ( mustanglm90_avg.dat' 
mm=assoc( 1 ,fltarr(512)) 

openr , 2 , ' [mustang ] m90_sdm . dat ' 
ss=assoc( 2, fltarr ( 512 ) ) 

openu , 3 ,'[ mustang ]m90  int.dat' 
dat=assoc ( 3 , f 1 ta r r ( 51I ) ) 

openu , 4 , ' [ mustang ]m90  intsdm.dat' 
dev=assoc( 4 , f Ita r r ( 51I ) ) 

openr , 5 ,'( mack ]mack  cal. dat' 
sen=assoc(  5  ,  f  1  ta  r  r  ( ■?!  2  )  ) 
sens=sen( 22 ) 
boolean= fltarr ( 512 ) 
dat ( 0 ) =mm( 0 ) 
dev ( 0 ) =ss ( 0 ) 

for  i=l,25  do  begin 

;get  array  of  where  dropped  words  are 
points=mm( i ) 

index=where ( points  eq  -1,0) 

; correct  for  sensitivity 

points=point s/sens 

, -assign  average  for  dropped  data  values 

points ( index )=( points! index+1 )+points( index- 1 ) ) /2 
, -correct  averages  for  double  drops 
boolean ( index ) =1 

double=where ( boolean* shi f t ( boolean, -1 ) ) 

points ( double )=( 2 *points ( double-1 ) +points ( double+2 ) )/3 . 

point s( double+1 )=( points (double-1 )+2*points(double+2 ) ) /3 . 

:assign  value  of  data  to  dropped  word  sdom's 
er rors=ss ( i ) 

errors! index ) =points ! index) 

:zero  first  six  elements 
points ! 0 : 5 ) =0 . 0 
errors! 0 : 5 )=0 . 0 
assign  data  to  m90_int.dat 
dat ! i ) =poi nts 

correct  for  sensitivity  and  assign  sdora  to  m90  intsdm.dat 
dev! i )=errors/sens  “ 

■ndf  or 


lose , 5 
lose , 4 
lose , 3 
lose , 2 
lose , 1 

eturn 

nd 


102 


;Walden ,1991 
pro  units 

openr , 1 ,'[ mustang ] m90_avg , dat '  ;reads  out  avg  spectra,  shifted 
a=assoc(3 ,fltarr(512)  ) 

openr , 4 ,'[ mustang ] m90_sdm . dat '  ;reads  out  SDOMs,  shifted 

c=assoc(4,fltarr(512)  ) 

openr , 2 ,'[ mack Imack  cal. dat'  ;calibration  to  correct  units  to  intensity 
b=assoc(2,fltarr(5ll) ) 
sens=b ( 22 ) 

openw, 3 , ' [ mustang ] m90_int . dat ' , 512* 4  ;file  to  accept  corrected  values 
int=assoc(3,fltarr(512) ) 

openw ,5, ' [mustanglm90_intsdm.dat' ,512*4 
intdom=assoc(5,fltarr(512) ) 

for  i=l,25  do  begin 

int ( i )  =  a ( i )/sens 
intdom( i )=c( i )/sens 

endf or 

close ,  1 
close ,  2 
close ,  3 
close ,  4 
close  ,  5 

re  turn 
?nd 
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; Walden , 1991 

;This  program  assigns  zeros  to  the  first  six  values  of  an  array  of  512 
pro  zerosix , fname 

openu , 1 , fname 
a=assoc(l,fltarr(5l2) ) 

temp=f Itar r ( 512 ) 

for  i=l,25  do  begin 
temp=a ( i ) 
temp( 0 : 5  )  =  0 . 0 
a ( i ) =temp 

endf or 

close  ,  1 

return 

end 
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APPENDIX  D 


Appendix  D  shows  the  programs  used  to  convert  the  scale 
factors  obtained  from  the  fitting  routine  to  absolute 
intensities.  The  equations  were  derived  from  Barth  (1965). 

The  equation  used  for  converting  the  scale  factor  to 
intensity  for  the  VK(0,5)  band  is: 

Iq  ^-VKCALCiO  .S)  *SF(al  t)  , 

where  Iq ^  is  the  band  intensity  in  Rayleighs,  VKCALC(0,5)  is 
the  total  intensity  of  the  synthetic  VK(0,5)  band  in 
Rayleighs,  and  SF(alt)  is  the  VK  scale  factor  for  the  altitude 
desired.  The  factor  VKCALC  is  just  the  area  under  the 
synthetic  spectrum  curve  for  the  transition  and  is  obtained 
using  the  IDL  command  TOTAL()  and  multiplying  by  wavelength 
bin  size,  as  seen  in  the  program  VK_ABSINT. 

The  conversion  for  the  LBH  scale  factor  included  all 
bands,  and  is  performed  in  the  program  LBH_ABSINT  and  the 
subroutine  INTENCALC.  The  basic  equation  for  this  program, 
derived  from  equations  (4)  &  (8)  of  Barth  (1965)  is  given  by: 

LBHCALCij ,  k)  Y, 

T  _ 

LBH  q  ' 

where  is  the  total  intensity  of  LBH  emissions  in 

Rayleighs,  LBHCALC(j,k)  is  the  area  under  the  synthetic 
spectrum  curve  for  a  transition  band  (j,k),  q  .  is  the  Franck- 
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Condon  factor  for  the  (j,k)  band,  and  v  is  the  corresponding 
frequency  for  each  transition.  The  summation  over  v"  is  to 
include  all  possible  upper  states.  The  result  of  this 
equation  was  essentially  the  same  for  any  (j,k).  This  thesis 
utilized  the  (4,12)  transition. 
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;Walden ,1991 

;This  program  converts  the  scale  factor  from  the  nimpfit_all  routine  into 
;an  absolute  intensity  for  a  particular  transition  band  of  Vega rd-Kaplan . 
;Stored  in  vk_absint.dat,  the  values  are  for  the  (0,5)  transition  in  record  1, 
;the  (0,6)  in  record  3,  and  the  (1,8)  in  record  5. 
pro  vk_absint 

~  openr , 10 , ' f inalpara .dat ' 
p=assoc(io,fltarr(25)  ) 
temp=p ( 0 ) 
alt  =  p( 21  ) 
vksf  =  p ( 1 3  ) 
delvksf =p( 14  ) 

openu , 11 , ' vk_ab5int .dat ' 
dat  =  assoc(ll,fltarr(25)  ) 
da t ( 0  )  =al t 

wl=indgen( 512 ) *3 .133+1798 
hold=f 1 ta  r  r ( 2  5 ) 
errhold=fltarr(25) 


for  run=2,6,+2  do  begin 

if  (run  eq  2)  then  begin  ;F-C  factors  are  not  currently  used 

vl=0  ;in  the  determination  of  intensity, 

v2=5  ;butareincludedforfutu reuse 

ql=5.9e-4  ;F-C  factor  for  (0,0) 

q2=1.69e-l  ;F-C  factor  for  (0,5) 

twk  =  2 . 6 

endi  f 

if  (run  eq  4)  then  begin 
vl  =  0 
v2  =  6 

ql=5.9e-4  ?F-C  factor  for  (0,0) 

q2=1.89e-l  ; F-C  factor  for  (0,6) 

twk=2 . 6 
endi  f 

if  (run  eq  6)  then  begin 
vl  =  l 
v2  =  8 

ql=3.32e-3  ;F-C  factor  for  (1,0) 

q2=7.87e-2  ; F-C  factor  for  (1,8) 

twk  =  l .  5 
endi  f 

for  a=0,24  do  begin 

int=indgen( 512)*0.0 

vksyn ,temp(a),vl,v2,wl,int 

int=int*twk 

hold(a)=total(int)*3.133*vksf(a) 

errhold(a)=total(int)*3.133*delvksf(a) 

endf or 

dat ( run-1 ) =hold 
dat( run)=errhold 

endf or 

1  ose  ,  1 0 
lose  ,  1 1 


e  tu  r  n 
nd 

[ synspec ] vksyn 

107 


, -Walden ,1991 
pro  lbh_absint 

;  This  program  determines  the  absolute  intensity  of  LBH  emissions  based 
;on  the  fits  of  intensity  to  three  different  transition  bands  (the  results  were 
, -basically  equal  in  Walden's  analysis).  Stored  in  record  1  of  lbh_absint.dat 
;is  based  on  th  (2,9),  record  3  is  the  (4,12),  record  5  is  (3,10). 
ope nr, 10, 'finalpara.dat' 
p=assoc(10,fltarr(25)  ) 
temp=p ( 0 ) 
lbhsf-p( 9 ) 
dellbhsf=p( 10 ) 
alt=p(21) 


openr , 11 , ' lbh_tots . dat ' 
t=assoc(il,fltarr(20)  ) 
lbhtotint  =  t ( 0  ) 

openu,12, 'lbh_absint.dat' 
i=assoc(i2,fltarr(25)  ) 
i  ( 0  )  =al t 

hold=f Itar r ( 25 ) 
errhold=fltarr(25) 

for  run=l,3  do  begin 

if  (run  eq  1)  then  begin 
v)  =2 
v2  =  9 
endi  f 

if  (run  eq  2)  then  begin 
vl  =  4 
v2  =  12 
endi  f 

if  (run  eq  3)  then  begin 
vl  =  3 
v2  =  10 
endi  f 

for  a=0,24  do  begin 

intencalc,temp(a) ,vl,v2,bandint 

inti bh=bandint*3. 133*1 bhsf(a)/lbhtotint( { temp(a)/100)-l ) 
errintlbh=bandint*3.133*dellbhsf(a)/lbhtotint((temp(a)/100)-l) 
hold ( a ) =intlbh 
errhold(a)=errintlbh 

endf o  r 

i ( 2* run-1 ) =hold  ;put  absolute  intensities  in  appropriate  record 

i ( 2  *  run ) =e r rhold  ;put  errors  in  correct  record 

■ndf  or 

lose ,  10 
lose ,  1 1 
lose , 1 2 

eturn 

nd 

[  synspec ] Ibhsyn 


108 


;  Walden ,1991 

;This  procedure  is  called  by  the  lbh_absint  program  in  order  to  determine 
; the  intensity  under  a  particular  "hump"  (a  transition  band)  in  the  LBH  curve. 
;The  inputs  vl  and  v2  are  the  upper  and  lover  states'  vibrational  levels 
;for  the  transition  to  be  determined.  Intlbh  returns  the  intensity  value, 
pro  i ntencalc , T , vl , v2 , intlbh 

openr,l, ' [synspecjfrnkcondn.dat' 
fc=assoc{i,fltarr(15,30) ) 

Ibhf c=f c ( 0 ) 

openr,2,'[walden] Ibht  rans_wl . dat ' 

.  wl  =  assoc ( 2 , f 1 ta r r  (  1 5  ,  3 0 ) ) 

1  bhwl=wl ( 0 ) 


summa=0 . 0 

wl=indgen( 512)*3. 133+1798 

int=fltarr(512) 

lbhfreq=fltarr(15,30) 

for  j=‘0,14  do  begin 

for  k=0,29  do  begin 

if  ( Ibhvl ( j , k )  ne  0)  then  Ibhf req ( j , k ) =2 . 997e8/lbhwl ( j , k ) 
if  (lbhwl(j,k)  eq  0)  then  Ibhf req ( j , k ) =0 . 0 
endf or 

endf  or 

for  i=0 , 29  do  begin 

summa  =  summa+ ( lbhfreq(vl, i ) *3)*lbhfc(vl ,  i  ) 

endf or 

lbhsyn,T,vl,v2,wl,int 
intvlv2=total( int) 

i n t lbh  =  i n tvl v2 * summa/(  Ibhfreq(vl,v2)''3*lbhfc(vl,0)*lbhfc(vl,v2)) 

close ,  1 
close ,  2 
return 
end 

synspecllbhsyn 
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divided  into  40oA  increments,  from  ISOOA  to  340oA,  and  are  in 
order  of  increasing  altitude.  Each  plot  shows  the  fit  versus 
data  in  the  upper  view,  and  a  breakdown  of  the  fit  by  major 
contributors  in  the  second  view. 


; Walden, 1991 

prc  plotapp,alt,xrain,XiT,ax 


openr,l,  'finalpara.dat:' 
pata  =  assoc(l,fitari-(25)  ) 

h9ight  =  para ( 21  ) 

openr,2,  '  [mustang]m90_intp.h.dat' 
d  =  assoc(2,fltarrl512M 

e  p e  n  r  ,  3  ,  '  f  i  n  a  1 1  wk  .  d  a  t ' 
f  =  assoc ( 3 , f Itar  n  512  )  ) 


tenip  =  pa  ra  (  0  ) 

it=( temp(alt-l ) )/100 

c  penr , 4 , 'allvk.dat' 
vk  =  assoc ( 4 , f 1 ta  r  r ( 512 )  ) 
pvk  =  pa  ra (13) 

0 penr, 5, 'all_lbh.dat' 
lbh=assoc(5,fltarr(512)) 
plbh  =  pa ra ( 9  ) 

cpenr , 6 , ' nodi t .  dat ' 
dlt  =  assoc( 6 , fltarr( 512  )  ) 

openr,7, 'noeps.dat' 
eps  =  assoc(7,fltarr(512)  ) 

openr,8, 'noqammack.dat' 
'jani  =  assoc  (  8  ,  f  1  ta  r  r  {  512)  ) 


pno  =  para( 11  ) 

epenr ,9 , 'OIII _:853 .dat' 
eiii=assoc(9,bltarr(512)  ) 

DO i i i =pa  r  a ( 7  ) 

:.ppnr,10,  'Oil_2470.dat' 
oii  =  assoc( 10,fltarr( 512)  ) 
poii=para(5) 

T penr, 11, '01  2972.dat' 
j  i  =  a  s  s  o  c ( 11, ? 1 1  a  r  r ( 512)  ) 
e  0  1  =  p  a  r  a  (  3  ) 

: p e  n  r  ,  1  2  ,  '  r;  1 1  ^  2  1  4  3  d  0 ij  b  .  d  3  t ' 

1 2  =  a  s  s  o  c ( 12, f 1 1  a  r  r i 512) ) 

;  n  2  =  p  a  r  a ( 15 ) 

1  =  i  n  d  g  e  n  \  512  '  *  3  .  1  3  3+  1  79  3  .  0  3  3 


;program  plots  data  vs.  fit  and  fit  by 
;parts  in  format  used  by  Walden  thesis 
! fc  r  appendix,  alt  can  be  1  to  25,  xmin 
;and  xmax  are  presumed  to  be  1800-2200, 
;  2200-  2600 , 2500-3000  ,  or  3000-3400,  so 
, -other  combinations  would  need  adjust- 
;ments  to  the  plot  key 


;data  is  in  records  d=l  to  25 


;fits  are  in  records  f=l  to  25 


i  I 


set  plot , ' hp ' 

device, /portrait ,/plotter_on_c  f f 
device, xoffset  =  l.  27 
device, yoffset=. 635 
device,xsize=17.78 
device, ysize  =  24. 13 

;TOP  OF  PAGE — DATA  AND  FIT 
set_xy , xmin, xmax 
set_viewport ,  .18,  .82,  .59,  .91 

if  (alt  ge  18)  then  $ 

!  rati  tle=  '  Data  &  Fit  at  ' +s  t  r  t  r  im  (  (  f  ix  (  he  i  ght  (  al  t-1  )  )  )  ,  2  )  +  '  km  (downleg)' 
if  (alt  It  18)  then  $ 

1 mti tle= ' Data  &  Fit  at  ' +s t r t r im ( ( f ix ( hei ght ( al t-1 ) ) ) , 2 ) + ' km  (upleg)' 

!xtitle='Wavelength(A) ' 

!ytitle='Intensity( R/A ) ' 

da  t  =  d (alt) 
f i t  =  f ( al t ) 

! 1 i ne  type  =  0 
plot , wl , da  t 

! linetype=2 
oplot , wl , f i t 

;KEY  FOR  TOP  PLOT--xrange  must  be  determined  in  term.s  of  pixels  so  that 
;yrange  can  be  known  for  placing  the  key  in  correct  place  on  the  graph 

;case  section  does  this  for 

case  1  of 

(xmin  eq  1800)  and  (xmax  eq 
ldat=dat (0:129) 
lfit=fit(0:129) 
end 

(xmin  eq  2200)  and  (xmax  eq 
ldat=dat( 129 ; 257  ) 

If it  =  f it( 129  :  257  ) 
end 

(xmin  eq  2600)  and  (xmax  eq 
ldat=dat( 257  :  384  ) 

If it  =  f it( 257  :  384  ) 
end 

(xmin  eq  3000)  and  (xmax  eq 
ldat=dat (384:511) 

If it  =  f it(  384  :  511  ) 
end 

else :  begin 

Ida  t  =  da  t 
1 f i t  =  f  i  t 
end 

endca  se 

:this  section  designates  x  and  y  coordinates  for  the  sample  lines  in  the  key 

xcoordl=[ xmin+20 ,xmin+50 ] 

xcoord2=[ xmin  +  20, xmin  +  100  ] 

xcoord3=[  xmin  +  120,xmiin+200  ] 

ycoordl=[  .96*max{ Idat)  ,  .96*max( Id  at)  ] 


the  ranges  der^o,:^*-8c 
2200):  begin 

2600 ) :  begin 

3000 ) :  begin 

3400):  begin 

;key  will  probably  not  be  correct 
;this  section  m.ust  be  manually  altered 
;to  concur  with  the  wl  range 
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ycoord2=[  .90*niax{ldat)  ,  .90*rr.ax(ldat)  ] 
ycoord3=[  .84*max(ldat) , .84*max(ldat)  ] 
yc oord4  =  [  .78*max(ldat) ,  .78*max(ldat)  ] 
ycoord5=[ .72*max(ldat), .72*max(ldat)] 
ycoord6  =  [  .66*max(ldat)  ,  .66*max(ldat) ] 
ycoord7=i .60*max(ldat) , .60*max(ldat) ] 

! 1 inetype=0 

•  oplot , xcoordl , ycoordl 
! 1 inetype=2 
oplot, xcoordl, ycoord2 
1 1 inetype  =  0 

xyouts,xmin+6C, ,995*ycoordl(0) , 'Data' 
xyouts,xmin+60, .995*ycoord2( 0) , ’ Fit' 

; BOTTOM  OF  PAGE — PLOT  OF  PARTS  AMD  FIT 

set_xy,xmin,xmax 
5et_viewport , .18, .82,  .09,  .41 

! mt i t le= ' Data  and  Fit  by  Parts' 

!xtitle='Wavelength{A) ' 

!ytitle='IntensityfR/'A)  ' 

!  1 inetype  =  0 
plot,v,-l,dat 

oplot, wl,dat  ;makes  data  a  bold  line 

!  1 inetype  =  l 

oplot, wl,vk(it)*pvk (alt-1)  ; dotted 

I  linetype  =  2 

oplot, wl,lbh(it)*plbh(alt-l)  ; dashed 

!  1 inetype  =  3 

nitric=pno(alt-l)*(gam( it)+eps( it)+.25*dlt(it) )/3.133e6 
,  oplot , wl , ni t r i c  , 'dash-dot 

! 1 i ne  type  =  4 

oplot, wl,oiii{0)*poiii(alt-l)  ; solid 

1 1 inetype  =  5 

-  oplot, wl,n2(l)*pn2(alt-l)  ; solid 

1 1 i ne  type  =  5 

oplot, wl,jii{0)*poii(alt-l)  ;solid 

! 1 i ne  type  =  5 

oplot,wl,oi(0)*poi(alt-l)  ;  solid 

;koy  for  bottom  plot--x  and  y  ranges  from  top  plot  used  for  locating  key 


:  a  s  e  1  of 

■xmin  eq  1800)  and  ( xmax  eq  2200):  begin 
! 1 i ne  type  =  0 
oplot ,xcoord2 , ycoordl 
oplot, xcoord2, ycoordl 
! 1 i ne  type  =  3 
oplot, xcoord2,ycoord2 
1 linetype=2 
oplot, xcoord2,y CO ord3 
! 1 i ne  type  =  5 
oplot,xcoord2,ycoord4 
! 1 1 ne  type  =  0 

xyouts,xmin-‘'110,  .99r*ycoordl^0)  ,  'Data' 
xyouts,xmin+110,  .995*ycoord2(G)  ,':;0' 
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xyouts,xmin+110, .995*ycQord3(0) , 'LBH' 

xyout s , xmin+1 1 0 , . 99 5*ycoord4 ( 0 ) , ' NI I  doublet  2143A' 

end 

(xmin  eq  2200)  and  (xmax  eq  2600):  begin 
! 1 i ne  type  =  0 
oplot,xcoord2,ycoordl 
oplot,xcoord2,ycoocdl 
! linetype=3 
oplot,xcoord2,ycoord2 
! 1 inetype=2 
oplot,xcoord2,ycoord3 
1 1 inetype  =  l 
oplot,xcoord2,ycoord4 
! 1 inetype=5 
oplot,xcoord2,ycoord5 
! i inetype=0 

xyouts,xmin+liO, .995*ycoordl(0) , 'Data' 
xyouts,xmin+110, .995*ycoord2(0) , 'NO' 
xyouts,xmin+110, .995*ycoord3(0) , 'LBH' 
xyouts,xmin+110, .995*ycoord4(0) , 'VK' 
xyouts,xinin+110,  .995*ycoord5(0)  ,  'Oil  247 OA' 
end 

(xmin  eq  2600)  and  (xmax  eq  3000):  begin 
1 1 i netype=0 
oplot,xcoord3,ycoordl 
oplot,xcoord3,ycoordl 
(line  type  =  3 
oplot,xcoord3,ycoord2 
! 1 inetype=l 
oplot,xcoord3,ycoord3 
! 1 inetype=5 
oplot,xcoord3,ycoord4 
! 1 i ne  type  =  4 
oplot,xcoord3,ycoord5 
! Iinetype=0 

xyout s ,xmin+210, .995*ycoordl(0) , 'Data' 
xyouts,xmin+210,  .995*ycooi:d2(0)  ,  'NO' 
xyouts,xmin+210, .995*ycoord3(0) , 'VK' 
xyouts,xmin  +  210, .995*ycoord4(0) , '01  291 2  A' 
xyouts,xmin+210, .995*ycoord5( 0) , 'OIII  2853A' 
end 

(xmin  eq  3000)  and  (xmax  eq  3400):  begin 
!  1 inetype  =  0 
oplot,xcoord2,ycoordl 
oplot,xcoord2,ycoordl 
!  1 inetype  =  3 
oplot,xcoord2,ycoord2 
!  1 i ne  type  =  l 
oplot,xcoord2,ycoord3 
!  1  i  ne  type  =  0 

xyouts,xmin+110, .995*ycoordl(0) , 'Data' 
xyouts,xmin+110, .995*ycoord2(0) , ' MO ' 
xyouts,xmin+110, .995*ycoord3(0) , 'VK' 
end 

’ 1 se :  dummy=  0 

?ndcase 


1  14 


;  c  1  o  5  e  ,  1  3 
r 1 ose , 1 2 
close , 11 


close  ,  10 
close  ,  9 
close  ,  8 
close  ,  7 
close  ,  6 
close ,  5 
cl ose ,  4 
close  ,  3 
c lose  ,  2 
close  ,  1 

return 

end 
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